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ABSTRACT 

We construct rest-frame luminosity functions at 3.6, 4.5, 5.8, 8 and 24/im over the 
redshift range 0<z<2 for galaxies and 0<z<4 for optical QSOs, using optical and 
infrared data from the Spitzer Wide-area InfraRed Extragalactic survey. The 3.6 
and 4.5/xm galaxy LFs show evidence for moderate positive luminosity evolution up 
to z ~ 1.5, consistent with the passive ageing of evolved stellar populations. Their 
comoving luminosity density was found to evolve passively, gradually increasing out 
to z~0.5-l but flattening, or even declining, at higher redshift. Conversely, the 24/im 
galaxy LF, which is more sensitive to obscured star formation and/or AGN activity, 
undergoes strong positive evolution, with the derived IR energy density and SFR 
density oc (1+z) 7 with 7=4.5tQ;g and the majority of this evolution occurring since 
z~l. Optical QSOs, however, show positive luminosity evolution in all bands, out to 
the highest redshifts (3<z<4). Modelling as L* oc(l+z) 7 gave 7=1.3^^ at 3.6um, 

7=1.0^01 a t 4.5u.m and stronger evolution at the longer wavelengths (5.8, 8 and 
24[im), of 7 ~ 3. Comparison of the galaxy LFs to predictions from a semi-analytic 
model based on CDM indicate that an IMF skewed towards higher mass star formation 
in bursts compared to locally is preferred. As a result the currently inferred massive 
star formation rates in distant sub-mm sources may require substantial downwards 
revision. 

Key words: galaxies: evolution - galaxies:photometry - quasars: general - cosmology: 
observations 



1 INTRODUCTION 

A landmark result in astronomy over the last two decades 
was the launch of COBE and the sub sequent discovery of 
a strong far-IR b ackground radiation llHauser et a.1.1 fl998l : 
iFixsen et al.lll99sl) . Studies of the extragalactic background 
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suggest at least half the luminous energy gener ated by stars 
has been reprocessed into the IR by dust (e.g. iPueet et alJ 
1996, Lag ache et al.lll999l) . suggesting that dust-obscured 
star formation was much more important at higher redshifts 
than locally. Space missions such as IRAS and ISO have 
resolved this background into large and significant popula- 
tions of IR sources - to luminous IR galaxies (Lir > 1O 11 L0) 
at z<l and to ultra-luminous IR galaxies (Lir > 1O 12 L0) 
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at z>l (e.g. see iDwek et alJ 119981 


Bareer et all Il999: 


iFlores et alJll99Sl ICharv & Elbazll2001 


: iFranceschini et alJ 


l200ll:lR,owan-R,obinson et alJl2004h. These strongly evolving 



populations reveal that dust-shrouded activity is pivotal in 
understanding the formation of stars and the central black 
holes in galaxies. The evolution that is required by such 
studies is in agreement with suggestions from longer wave- 
lengths - the luminous sub- mm galaxies have been shown to 
have median redshifts > 2 ^Hughes et al jflQQSl : IScott et alJ 
120021 : Ichapman et al.ll2005l) . indicating that significant evo- 
lution has occurred. 

Complementary results come from Hubble Space 
Telescope observations, particularly when combined with 
Spitzer observations. These show that the bulk of galaxy 
morpholo gies we see locally were put into place over 
0.5<z<3 JConselice et~ai]l2004l ; Ipapovich et al]l2005h . 

A combination of these results suggests that, in order 
to understand the mass assembly history of the stars and 
central black holes in galaxies, studies need to include IR 
observations and be made out to high redshift. 

A basic tool in the study of galaxy populations is the 
construction of Luminosity Functions (LFs). These have 
long been used to constrain galaxy formation models and 
to quantify star formation rates (SFR) and evolution (both 
of luminosity and number density) but have historically been 
focused on optical wavelengths. Given the importance of IR- 
based observations to galaxy evolution we, potentially, can 
make large strides in understanding how galaxies assemble 
by constructing near, mid and far-IR galaxy LFs, and seeing 
how they evolve with redshift. The shorter wavelengths tell 
us how evolved stellar mass builds up, whilst longer wave- 
lengths indicate when and where the observed starbursts 
that are forming these stars are occurring. 

To date this has been difficult since large, deep sam- 
ples are required in order to trace the evolution - earlier 
work using IRAS, ISO or ground-based sub-mm observa- 
tions have been hampered by sensitivity, area and (for sub- 
mm) identification difficulties. One overarching feature of 
these earlier studies has, however, been the strong evolu- 
tion demonstrated by their LFs. Some of the earliest IR 
LFs were constructed f r om th e IRAS satellite mission (e.g 
iRowan-Robinson et a l derived 25, 60 and 100u.m 
LFs~). ISaTmdljnH^tH!dj " ifl990l) used IRAS data to derive the 
60\im and 40-120(i.m LFs, which were found to be indica- 
tive of strong evolution such t hat luminosity increase s with 
redshift, perhaps cx(l+z) 3±1 . IClements et alJ (|200l|) used 
deep 12u.m ISO data with follow-up optical imaging and 
spectroscopy to determine the mjy 12\im LF, finding that 
the excess (in comparison to the low-redshift IRAS re- 
sults) at high luminosites (higher redshift) was compatible 
with rapid lum i nosity evolution oc(l+z) 4,5 . More recently, 
ISerieant et alJ i2004l) used the o ptical-IR band-merged 
ELAI S Final Analysis catalogue of iRowan- Robinson et all 
<2004h to calculate the ELAIS 90um LF, finding that, 
for consistency with source counts, a luminosity evolu- 
tion of (l+z) 3,4±1 was required - consistent with the evo- 
lution in comoving volume averaged SFR at z< 1 de- 
rived from rest-frame o ptical and UV surveys (e.g SLOAN; 
iGlazebrook et alJl2003h . 

The launch of Spitzer JWerner et alJl2004 has given 
a powerful tool to astronomers. Its wide-area deep-mapping 
survey capabilities means extensive progress can be made 



in our understanding of galaxy assembly, allowing the con- 
struction of LFs in multiple IR bands across wide luminos- 
ity and redshift range s. The Spitzer Wide-area InfraRed 
Extra galactic (SWIRE: lLonsdale et all2003l : lLonsdale et alJ 
120041) s urvey makes use of all 7 of Spitzer's imaging bands 
(IRAC. iFazio et al] I200I 3.6, 4.5, 5.8 and 8^m. MIPS, 
iRieke et all 12004 24, 70 and 160u.m) and extensive ancil- 
lary data at other wavelengths across a total area of ~ 49 
deg 2 and has detected around 2.5 million infrared sources. 
At the shorter IRAC wavelengths the survey can provide 
constraints on the more evolved stellar populations whilst at 
longer wavelengths it becomes increasingly sensitive to on- 
going star formation. As such, SWIRE is able to investigate 
the evolution of the stellar mass function, star formation 
rates and probe galaxy luminosity and density evolution. 
With the combination of large area and high sensitivity it 
is uniquely able to provide good statistics for the character- 
istic IR sources out to z>l whilst detecting rare, luminous 
sources out to high redshift. 

In this paper the power of SWIRE, combining both IR 
and optical data, is used to investigate galaxy and Type 
1 active galactic nuclei (AGN) luminosity functions so that 
questions on their formation, evolution and relationship over 
a wide redshift range can be addressed. 

To carry this out, the template-fitting photometric red- 
shift method is first extended to incorporate near-IR data 
and then applied to band-merged 3.6 to 24[im Spitzer data 
and associated optical photometry (U, g', r', i' and Z) from 
one SWIRE field (ELAIS-N1). The resulting photometric 
redshifts and best-fitting optical and IR template solutions 
are then used to obtain luminosity functions (LF's) and to 
investigate evolution with redshift. 

In Sj5] the IR and optical catalogues, data calibration 
and band-merging are set out. Photometric redshifts and 
luminosity function calculations are outlined in J3]and 
Luminosity functions with redshift are presented at 3.6, 4.5, 
5.8, 8 and 24(im for galaxies and Type 1 AGN in J5| and 
comparisons to other work are made in The results of 
fitting functional forms to these LFs, evolution with redshift, 
energy density and star formation history are reported in J7| 
and overall discussions and conclusions are presented in jJS] 
For this work the flat cosmological model with Ho=72km 
s~ 1 Mpc~ 1 and A=0.70 is used. Throughout this paper, the 
term "AGN" is used to refer to all Type I active galactic 
nuclei as determined from optical to near-IR template-fits. 



2 THE DATA 
2.1 Infrared data 

In this paper we consider one field observed in the first 
stage of the SWIRE survey. This is 6.5 deg 2 of the ELAIS- 
Nl field (specifically, the area in common with the Isaac 
Newton Telescope Wide-Field Survey (INT WFS) survey, 
see H2.21 . centred on 16/illm + 54d55m, for which we have 
band-merged 3.6, 4.5, 5.8, 8 and 24u.m catalogues. Fur- 
ther information on the IRAC and MIPS data processing 
can be found in t he SWIRE Data Release 2 Document 
jSurace et alJl2005h and in Shupe et al. (2006:in prep.). A 
plot of the IRAC, MIPS and optical coverage of ELAIS-N1 
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is shown in Fig. The depths (5 a) are approximately 3.7, 
5.4, 48 and 37.8 fiJy in the IRAC bands (3.6 to 8\im) and 
around 230^iJy at 2A\im (not e this is a factor ~2 higher than 
original estimates for 24u.m: ISurace et alj|2005h . 



2.1.1 Reliability and completeness 
Reliability 

One of the most important sources of unreliability are 
bright star artefacts. The majority of artefacts associated 
with bright stars were eliminated from the basic SWIRE 
source lists by tuning SExtractor to avoid their close 
vicinity using bespoke masks. Visual inspection of random 
subsets of the catalogue were carried out at all wavelengths, 
with results indicating a reliability > 99%. Investigation 
of the observed SEDs of SWIRE sources as compared to 
those of known o bjects was also carr ied out as a second 
test of reliability feurace et alJ 12005). In order to assess 
the reliability of the MIPS 24u.m detections, Shupe et 
al. (2006:in prep.) have compared the SWIRE data with 
deeper GTO data processed using the same methods and 
matching extracted sources. Defining reliability as the ratio 
of all matched sources to all extracted SWIRE sources 
gives a reliability better than 98% for fluxes above 0.35 mjy. 

Based o n studies with the band-merge d optical-IR cat- 
alogues (e.g. [R owan- Robinson et alj|2005l hereafter RR05; 
ISurace et a n l2005h . a reliable source is here considered to 
be one that has a detection with S:N>5 in the band in 
question and a further detection in another optical or IR 
band. Although this approach will conceivably exclude real 
sources which are only detected in a single band this is not 
likely to be a significant population. The largest effect might 
be expected for the 2A\im population - if there is a pop- 
ulation distinct from the shorter wavelength sources then 
the requirement of detection in another band will introduce 
a bias. Of the 40,387 24u.m sources with S:N>5, 772 have 
no detection at any other wavelength in the catalogue, or 
<~^2 per cent. Thus any introduced bias is not severe (note 
that a more likely source of bias will be introduced later 
by requiring that the source obtains a photometric redshift). 

In the area in common with the optical data (sec H2.2H 
331,209 Spitzer sources (81%) are defined as reliable. Num- 
ber counts are shown for each band in Fig. [5] 

Completeness 

The completeness of the IRAC data was evaluated by 
ISurace et all (120051) using a de eper field wit hin ELAIS-N1 - 
the GOODS Validation Field (ICharv et alJl2004l) . taken as 
part of the Extragalactic First Look Survey (FLS) program. 
FLS and SWIRE data were compared in a 25 square ar- 
cminute area. The SWIRE 95% completeness level was cal- 
culated to be at 14/xJy, 15/xJy, 42/xJy, and 56/iJy for IRAC 
3.6u.m 4.5u.m 5.8u.m and 8[im respectively. A secondary 
analysis via differential source counts agreed well. The over- 
lap of 2A\im SWIRE and FLS sources was not large enough 
for meaningful statisticsso in order to derive the MIPS 2A\im 
completeness function, simulated sources were added to the 
central high-coverage portion of the ELAIS-N1 mosaics and 



SExtracted, allowing completeness versus flux to be calcu- 
lated (for more details see Shupe et al. 2006: in prep). The 
completeness level at 500^iJy was calculated to be ~97%. 
Shupe et al. also carry out a direct estimate by comparing 
SWIRE data to deeper GTO data processed in the same 
manner, finding good agreement with the simulations for 
the drop in completeness with decreasing flux. 

2.1.2 Bandmerging 

Bandmerging of IRAC and MIPS 24\im catalogues was car- 
ried out with the Spitzer Science Center's bandmerge soft- 
ware. Bandmerge reliability was es timated from cumu lative 
distributions of positional offsets JSurace et alJ HOPS'). For 
> 90% of merges, sources within 1.5" for IRAC pairs or 
within 3" for IRAC-MIPS-24 can be considered reliable. 

2.2 Optical data 

In this paper we use optical photometry available in the 
ELAIS-N1 field (U, g', r', i' and Z) from the INT WFS. 
For more information on the survey and its goals see 
iMcMahon et alJ i200ll) . Basic reduction, photometry and 
catalogue information is outlined below: 

The ELAIS-N1 (16hl0m+54d30m) field is approxi- 
mately 9 deg 2 in U, g', r', i' and Z bands. These data have 
been processed at Cambridge Astronomical Survey Unit's 
reduct ion pipeline (CASU) as described in llrwin fc Lewisl 
(2001). 

The mean limiting magnitudes (Vega, 5a in 600 sec- 
onds) in each filter across this area are 23.40(U), 24.94(g'), 
24.04(r'), 23.18(i'), and 21.90(Z), however the completeness 
falls off before this, as is demonstrated by the r' magnitude 
distribution of sources in Fig. [5] 

2.2.1 Photometric calibration and merging 

Data reduction was carried out with the CA SU pipeline (fur- 
ther details are available in iBabbedeeteOol and S urace et alJ 
2005) with an overall photometric calibration for the survey 
at the level of 2 per cent. In order to merge the individual 
optical pointings into a single optical catalogue an associa- 
tion radius of 0.9" was used (based on the seeing, 1.17", 
and a nearest-neighbour separation histogram). For sources 
that had duplicates (from overlapping tiles), the one with 
the lowest photometric error was retained. 

2.3 Optical associations 

The Spitzer and optical observations do not have a 100% 
overlap. For ELAIS-N1, the area overlap is ~ 70%, as 
illustrated in Fig. The optical catalogues were asso- 
ciated with the band-merged SWIRE 3.6 to 2A\im cat- 
alogues using the Webcmp package of IPAC's Infrared 
Science Archive (IRS A). Webcmp is optimized for very 
fast cross-correlations, by position, between extremely large 
datasets. A radius of 1.5" was used for cross-identifications, 
a value chosen to maintain completeness whilst keeping false 
matches to a minimum. Where an IR source picks up mul- 
tiple optical associations, the association with the highest 
probability indicator value, an internal Webcmp statistic 
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Figure 1. The Spitzer IRAC data (grey) and the overlap area 
with optical data (black). The 'crenellations' are due to the layout 
of the WFC CCDs. The MIPS observations are indicated by the 
grey dotted region, almost entirely obscured by the IRAC area. 
Image from Surace et al. 2005. 



based on the match algorithm of the NASA/IPAC Extra- 
galactic Database (NED), was retained. This match statistic 
is, to first order, based on the angular separation, weighted 
by the catalogue uncertainties. 

Completeness and reliability of these S pfeer-optical 
cross- identifications (XIDs) was investigated in lSurace et alJ 
( 2005) which showed that out to the maximum association 
radius of 1.5" the Spitzer-optical XIDs are essentially 100% 
complete, excepting the possibility of a few bright true 
matches at large separations. The ratio of 'false' matches 
to 'real' matches suggest the reliability is ~ 95%. 

Of the 331,209 reliable Spitzer sources in the overlap 
area, 82,606 are optically blank (including a small fraction 
in INT WFS coverage gaps), meaning 248,603 (75%) of the 
reliable Spitzer sources have optical counterparts. We then 
make the further requirement, in order to obtain reliable 
photometric redshifts in H3.1.3I that the source be detected 
in four of the bands U, g', r', i', Z, 3.6[im, 4.5u.m which 
reduces the sample to 180,879 sources (74%). 

2.4 Star/galaxy separation 

Star/galaxy separation is an important but complex issue, 
even with the luxury of multi-wavelength data. Stellarity 
in a band gives a measure of the compactness of a source, 
however this measure is not an infallible star identifier - at 
faint fluxes (up to 2.5 magnitudes above the limiting mag- 
nitude) the measure begins to break down since one can 
then no longer reliably differentiate between extended and 
point-like sources. Even at brighter fluxes, a star may not 
be classed as point-like in all bands (e.g. saturation effects) 
whilst conversely, though a nearby galaxy will be clearly ex- 
tended, distant galaxies and quasars will also appear point- 
like. Thus effective discrimination requires consideration of 
stellarity and saturation flags, magnitudes and colours. 

For the identification of stars at fainter magnitudes it is 
worthwhile to note that at wavelengths covered by IRAC 
we are well into the Rayleigh- Jeans tail of stellar black- 
body curves and that stars all have essentially the same 
slope in their SEDs. This means that they will be colour- 



less in the IRAC bands (in Vega magnitudes). As such, the 
stellar locus is well-defined in colour-colour space and can 
be used in conjunction with stellar flags to define a robust 
star/galaxy separator. Here, the (3.6u.m/r') versus (r'-i') di- 
agram in conjunction with stellarity is used as the basic 
method to identify stars and remove them from the cata- 
logue, as detailed in RR05. From these definitions, 41,042 
of the 180,879 optical-IR sources are classed as stars (23%), 
leaving 139,823 sources. 

For the optical-IR catalogue, almost all sources with 
r' < 17 are stars and hence for the analyses in this paper we 
assume that all sources brighter than this are stars and re- 
moved. In addition, all sources that are flagged as saturated 
in one of more optical bands are removed, since in addition 
to almost certainly being stars, their photometry and stel- 
larity designations will be unreliable. As a final cut, those 
sources fainter than r'=23 are removed in order to main- 
tain reliability and completeness (see Fig. 0. Thus we are 
left with a catalogue of 102,645 extragalactic sources, with 
17 <r' < 23. These sources can now be split into extended 
and point-like sources, where point-like sources will be either 
quasars or faint galaxies. Sources are defined as point-like 
(and so Type 1 AGN templates are considered in addition 
to the galaxy SEDs in photometric redshift fitting as will be 
described in H3.1H based on the optical cla ss flags. The best 
bands, as found in iBabbedee et alJ J2004I) (hereafter B04), 
to define stellarity are g', r' and i' so if the flag is set to 
'point-like' in any of these bands then the source is defined 
as point-like. The goal here is simply to exclude those sources 
that do not look point-like in any band from having AGN 
templates applied since this cuts the number of spurious 
high-redshift AGN fits to low-redshift galaxies. The r' mag- 
nitude distribution of sources defined as stars is compared 
to that of the remaining sources in Fig. |2] and it can be seen 
that the star distribution is relatively flat with magnitude, 
and dominates at the bright end. Th e stellar counts from 
the stellar population synthesis code of lGirardi et al] (120051) 
for their CDFS field have been overlaid, converting from R 
band to r' and scaled to the same area. Both the CDFS 
and our ELAIS-N1 field are 'extragalactic' fields, away from 
the Galactic plane and bulge, so have similar proportions of 
each stellar populatio n - a normal Galac tic stellar population 
sequence fELAIS-Nl: llbata et aT1l2003() . We apply an addi- 
tional renormalisation factor (0.74) to the total counts in the 
Girardi et al. model to mimick the cut placed on our multi- 
band catalogue prior to the identification of stars. This cut 
was the requirement of detection in four of the bands U, g', 
r', i', Z, 3.6u m, 4.5u.m, which reduced our sample by a factor 
0.74. The assumption that this same factor can be applied 
to the Girardi et al. counts appears reasonable since the Gi- 
rardi et al. counts then match our observed stellar counts at 
~r'=17, where we can be confident our stellar classification 
is accurate. The open squares mark the approximate point at 
which incompleteness starts to hav e an effect at the faint end 
(based on the observed counts of iGroeneweeen et al ] |2002| 
in the CDFS). This empirical comparison shows that the 
observed star counts distribution is in agreement with that 
expected from an extragalactic field. At the faint end where 
our ability to identify stars is less complete, some may be 
passed to the redshift code. However, the majority of these 
appear to end up in the histogram of sources that fail to 
obtain a redshift - comparison of the summed distribution 
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Figure 2. Number count distributions. First panel: Magnitude distribution (Vega) of r' band sources in the ELAIS-N1 full-coverage 
region (full optical-IRAC-MIPS overlap). Black histogram - all sources; cyan - those that obtained a photometric redshift solution; green 

- those that failed to obtain a redshift solution; red - those sources defined as stars and removed; purple - the sum of stars and photometric 
redshift failures. Vertical dashed lines indicate the nominal 5<r limit of 24.04 and the magnitude cuts at r'=23 and 17. The filled circles 
joined by a dotted line represent the stellar counts from the stellar population synthesis code of Girardi et al. (2005), empty squares are 
where incompleteness will have a dominant effect. Remaining panels: Number counts in the full-coverage region for IRAC and MIPS 
24um sources. Black histogram - all reliable sources with a detection (S:N>5) in that Spitzer band (will include stars); red histogram 

- the subset that are optically blank (will include stars); green histogram - the subset that have optical detections (will include stars); 
cyan histogram - those with an r' band detection of 17< r' < 23 and which obtained a photometric redshift with X^ed < 10- Vertical 
red dash- triple-dot lines - chosen flux below which to cut sources from the LF calculations in order to maximise completetness. 



of stars and photometric redshift failures accounts for the 
majority of 'missing' stars at the faint-end. It is concluded 
that the contamination of stars incorrectly fit as galaxies or 
AGN at some redshift is small. 



3 TEMPLATE FITTING - A TWO-STAGE 
APPROACH 

The emission of a galaxy at optical to near-IR and at mid-IR 
wavelengths is not necessarily well correlated - that is, the 



optical to near-IR emission will often be dominated by emis- 
sion from stars, whereas at longer wavelengths re-processed 
emission from dust or synchrotron emission may be the dom- 
inant contribution. Hence for a single emission mechanism 
the optical to mid-IR SED of a galaxy cannot be allocated 
to one of a set of generic templates. Here, the analysis is 
separated into optical/near-IR and mid-IR, and template- 
fitting is carried out in a two-stage approach as in RR05: 
Optical and IRAC 3.6utti, 4.5\im detections are used for 
the photometric redshift fitting procedure (see Section foMH : 
the longer wavelengths are then used to fit IR templates to 
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those sources defined as having an IR-excess (see Section 

3.1 Photometric redshifts and extension to IR 
bands 

In order to calculate LFs the redshifts of the sources are re- 
quired. Here the updated version of the ImpZ code of B04 
is extended to utilise some of the IR bands from Spitzer 
and demonstrated on several sets of spectroscopic redshifts 
f A3. 1.21 . We seek to demonstrate both the validity and im- 
provement gained by extending the template fitting photo- 
metric redshift technique to incorporate near-IR data, prior 
to applying it to the full ELAIS-N1 optical- Spitzer cata- 
logue. A more in-depth study of the photometric redshift 
technique, its extension to IR and application to a number 
of large spectroscopic samples will be presented in future 
work (Rowan-Robinson et al. 2006 in prep.; Andreon et al. 
2006 in prep.). 

An important requirement for photometric redshift fit- 
ting is that the photometry in each band is measured in the 
same way - different aperture sizes will sample the galaxy 
out to different radii, and possibly in different locations. 
When only using optical data, use of fixed aperture fluxes 
is sufficient, however if we include IRAC near-IR fluxes it is 
important to be comparing like with like, and these fluxes 
are integrated fluxes. Hence here, for the optical bands, the 
aperture magnitudes corrected to the integrated magnitudes 
via curve-of-growth analysis are used. This procedure could 
introduce errors for galaxies whose integrated SEDs differ 
largely from the SEDs of their central regions, but in prac- 
tice the results are consistent, as found in A3. 1.21 

3.1.1 Photometric redshift method 

Improvements to the code include an iteratively calculated 
correction to the calibration in each band (suggested by 
lllbert et alJl2006l) and a parabolic interpolation of the red- 
shift solution space to improve the precision, along with im- 
provements based on studies in RR05. One issue raised by 
RR05 was that for some AGN there is a significant dust 
torus contribution to emission at 3.6 and 4.5u.m which is 
problematic if the code attempts to fit it as stellar emission. 
One important change to the ImpZ implementation to the 
SWIRE data is to carry out a double pass of the code on the 
data, in order to deal with AGN dust tori successfully: In 
the first pass, the 3.6 and 4.5u.m data are not included in the 
fit if S(3.6)/S(r') >3; in the second pass the 3.6 and 4.5u.m 
data are included in the fit provided S(3.6)/S(r') <300, ex- 
cept when the mid-IR fitting resulted in an AGN dust torus 
fit. Mid-IR excess will instead be fit separately by more ap- 
propriate SEDs - see H3.2I 

The chosen ImpZ code setup parameters were as follows: 
Templates: The six galaxy templates and two (Type 1) AGN 
templates of B04 are shown in Fig. [3] The two AGN tem- 
plates only differ from one another at ~ 3000A where one 
has a 'small blue b ump' in its emissi on (blend of Balmer con- 
tinuum and Fe IL lMaoz et al.llT993r) . Determining the exact 
shape of the restframe far-UV spectra of AGN is problem- 
atic. Here a simple continuation of the power-law is adopted, 
though at redshifts where this region enters the optical fil- 
ters, the dominant effect is due to IGM absorption. IGM 
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Figure 3. The templates used by ImpZ, showing the galaxy 
templates in green (E, Sab, Sbc, Scd, Sdm and Sb galaxies) and 
the two AGN templates in purple. 



treatment jMadau et al] ll996l and Galactic extinction cor- 
rections are included by ImpZ. 

Internal extinction: Variable extinction is allowed in the fits 
via variable A„ added to the templates. Ay limits of to 
2 in the Ay freedom are used. As found in B04, the in- 
clusion of A v as a parameter gives an important improve- 
ment in the accuracy o f the resulting redshifts, though from 
iBabbedee et alj <l2005l) it is also apparent that the derived 
extinction of the source is imprecise. With the inclusion of 
near-IR data in the fits, A„ freedom is now extended to the 
AGN fits, though restricted to A„ ^ 0.3 in order to prevent 
aliasing with normal galaxies (it was found that otherwise 
low redshift galaxies could be fit as higher redshift AGN 
with high extinction). 

Ay Prior: A prior expectation that the probability of a given 
value of Ay being 'correct' declines as |Ay| moves away 
from is applied. This is introduced by minimising Xred 
+ a Ay rather than \ 2 (tt=3 here). This use of a prior can 
be viewed as a weak implementation of Bayesian methods 
and was reached based on degeneracy studies in B04 and 
R owan- Robinson! (12003ft . hereafter RR03. 
Magnitude limits: RR03 and B04 found it necessary to ap- 
ply absolute magnitude limits to exclude unlikely solutions 
(such as super-luminous sources at high redshift). Here, the 
following limits are used: Absolute magnitude limits of - 
16<M B <[-22.5+0.3z] for galaxies and -21.7<M B <-26.7 for 
AGN. 

A 'good' photometric redshift solution is chosen to be one 
with a Xred < 10- This value of Xred is chosen as a threshold 
based on the results in B04 and this paper, since fits above 
this value are less reliable. 
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3.1.2 Comparison of spectroscopic and photometric 
redshifts 

Investigations have been carried out for several datasets in 
the SWIRE fields, all of which have optical photometry, IR 
data from Spitzer and spectroscopic redshifts. RR05 pre- 
viously demonstrated the success of the code for both an 
ELAIS-N1 sample and for a sample in the Lockman Vali- 
dation Field (the WIYN survey; Owen et al., in prep.). It 
is worth noting that the results from the WIYN compari- 
son showed the code is successful for galaxies out to z>l, 
whilst earlier work on HDF-N (RR03) and CDF-N (B04) 
samples have shown success for galaxies in the range 0<z<5 
and 0<z<2 respectively, and for Type 1 AGN out to z~5 
(B04). 

Here, updated results for several spectroscopic samples 
are outlined below, measuring the reliability and accuracy of 
the photometric redshifts via the fractional error Az/(1 + z) 
for each source and examining the mean error Az/(1 + z) 
and the rms scatter cr z : 

(i) LOCKMAN VF: In the 0.3 deg 2 LOCKMAN Vali- 
dation Field (VF) of the SWIRE survey there are U, g', r' 
and i' data, taken with the MOSAIC I camera on the May- 
all 4m telescope at Kitt Peak National Observatory, and 
Spitzer 3.6 to 24/im data. The comparison of photometric 
redshifts to 269 sources with spectroscopic redshifts from 
the WIYN survey (Owen et al., in prep.) gives an rms scat- 
ter of 0.061, with galaxies out to z~l, and Type 1 AGN 
out to z~3. Compare this to typical accuracies of o~ r ed ~0.1 
for studies considering only optical photometry in the fits 
(RR03). 

(ii) ELAIS- N1: The spectroscopic r edshift sample used 
in RR05 and iRowan- Robinson et alJ i2004h . comprising 
spectroscopic redshifts in ELAIS-N1 (reported by Perez- 
Fourn on et al. 2006: in prep, and Serjeant et al. 2006: in 
prep.; lHatziminaoglou et alJl2005bT) with optical U, g', r', 
i', Z and J, H, K (where available) and Spitzer data. For 
this sample the total rms scatter, a to t, was 0.057, where the 
majority of sources out to redshift z~0.5 are galaxies whilst 
at higher redshifts (to z~3) Type 1 AGN dominate. 

(iii) VVDS: Comprises optical photometry from the 
CFH12K-VIRMOS survey dLe Fevre et al.1 120041) : a deep 
BVRI imaging survey conducted with the CFH-12K camera 
in four fields. Additionally, there are U band data from the 
2.2 m telescope in La Silla and J, K data from the NTT. The 
spectroscopic redshi fts come from the fo llow-up VIRMOS- 
VLT Deep Survey dLe Fevre et al]l2003lL a large spectro- 
scopic follow-up survey which overlaps with the SWIRE 
data in the XMM-LSS field. For 1394 galaxies in the range 
0<z<1.5 the achieved accuracy was cr z =0.052 with very low 
systematic offset. 

For all the samples, the mean systematic offset between the 
photometric and spectroscopic redshifts was found to be es- 
sentially zero to the precision of the photometric redshifts. 
For example, the ELAIS-N1 sample has ~Kz/(l+z) = 0.0037. 
In comparison to earlier results for the spectroscopic samples 
where photometric redshifts were derived from only the opt i- 
cal photometry (RR03; B04: iRowan-Robinson et alJl20ol) . 
the inclusion of the 3.6 and 4.5u.m data from Spitzer en- 
ables several solution degeneracies to be broken, since data 
is provided on the underlying SED across a much wider ob- 



served wavelength range (0.3u.m to 5u.m). The additional 
near-IR information has reduced the dispersion and have en- 
abled the rejection of most of the extreme outliers resulting 
from optical-only results. This is an important improvement 
- even though the two IRAC bands may not have enabled 
a correct photometric redshift to be found for some of the 
original extreme outliers they have enabled the code to re- 
ject most of the original optical-only bad outliers, leaving 
fewer spurious redshift solutions. 

From these studies it is concluded that the inclusion 
of the first two IRAC bands in the photometric solution 
reduces the rms error and enables the rejection of the ma- 
jority of catastrophic outliers. Furthermore, the inclusion of 
the first two IRAC bands in the solution reduces the rms 
scatter. Further information on photometric redshifts and 
their extension to IR wav e length s for the SWIRE survey 
can be found in Babbcdgc (2004), RR05, Rowan- Robinson 
et al. (2006 in prep.) and Andreon et al. (2006 in prep.). 

3.1.3 Application to the optical- Spitzer catalogue 

The photometric redshifts of the optical-Spitzer 
extragalactic sources have been derived from their 
optical and IRAC 3.6, 4.5u.m detections using the updated 
version of ImpZ with the setup set out in H3.1.1I The 
resulting photometric redshifts are expected to be accurate 
to ~ 0.05 in logio(l + z), based on comparison to the 
spectroscopic sample in i)3. 1.21 and the work of RR05, B04 
and RR03. 

Due to the aliasi ng and photometric var iability issues 
noted in B04 (and also lAfonso-Luis et al.l2004f) the accuracy 
for the AGN will be lower, ~ 0.2 in Zo.<?io(l + z). Th e study 
of SWIRE/C/iandra sources bv iFranceschini et al.l (2005) 
(F05) also found that the p hotometric redshift ac curacy 
of their AGN was low, whilst IKitsionas et al.l (|2005) found 
that, for their X-ray selected XMM-iVeu>ton/2dF sample, 
the photometric redshift accuracy for AGN was < 0.2 in 
Az/(1 + z) for 75 per cent of the sample, with indications 
of the need for a new set of templates in order to model the 
outliers. Here, in order to successfully deal with cases where 
there is significant dust torus contribution to emission at 3.6 
and 4.5u.m, a 'double pass' of the catalogue through ImpZ 
is carried out as described in H3.1.1I 

Redshift results 

Of the 139,823 extragalactic sources input to the code 96% 
of sources obtain a photometric redshift solution. For the 
102,645 sources with 17 <r' < 23, there are only 800 red- 
shift failures and there are 96,177 with solutions below the 
Xred maximum of 10. It is this 17 <r' < 23 sample with 
redshift solutions of Xred < 10 that we will make use of in 
the calculation of LFs (see Tabled for statistics). 

From the plot comparing of the r' magnitude distribu- 
tion of sources that did and did not obtain a redshift (first 
panel of Fig.|5J one feature is apparent: the magnitude dis- 
tribution of those sources that failed to obtain a redshift 
is bimodal - there is a small distribution of bright sources 
(r' <20) which most likely failed to obtain a solution due to 
problems with saturated photometry and miss-associations 
of artefacts around bright stars and there is also a distri- 
bution of fainter sources (r' >21) for whom failure can be 
attributed to the lower accuracy of photometry at fainter 
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fluxes, the loss of bands towards the limiting depths and 
also to them being faint stars which were not picked out by 
the original star/galaxy separation. 

We plot Xred distributions of the photometric redshift 
solutions in Fig. [I] split into galaxy and Type 1 AGN tem- 
plate solutions. In each case, the distribution has been split 
into the same redshift bins as will be used for the various LF 
determinations in f0J and the typical proportion of galaxy 
and AGN solutions in each redshift bin with Xred < 10 i s 
> 90%. For the galaxy solutions the slope of the Xred distri- 
butions does not change with redshift, indicating that the 
applicability of the method and chosen templates remains 
successful out to the maximum redshift considered for the 
galaxy LF calculations. The slope of the Xred distribution 
for the AGN solutions does show some variation, with indi- 
cations that the solutions improve at the highest redshifts 
considered (3 to 4). This is attributed to the shifting into 
the optical bands of the Lyman break and IGM features in 
the AGN template which allows the AGN to be identified 
more readily. 

Galaxy solutions 
We plot redshift distributions of selected IR populations 
(3.6u.m and 24(j.m sources 1 ) in the left panels of Fig. 
El The distributions for 4.5 and 5.8p.m sources are sim- 
ilar to the 3.6[im whilst the 8u.m distribution starts to 
move towards the distribution seen for the 24\im sources. 
The galaxy populations, with (z ga ;)=0.55, extend out to 
a redshift of around 1 before dropping away. For compar- 
ison, the 24|xm-selected (to 83/xJy) galaxy population of 
iPerez-Gonzalez et alJ i2005l) (PG05) peaked at a redshift of 
of 0.6-1, and decayed monotonically from redshift 1 to 3. A 
c omparison of the galaxy r e dshift distributio ns t o the models 
oflRowan-Robinsonl feOOlft . lXu et alJ fcOOSt and lPozzi et"aH 
(2004) is given in RR05 and indicates that the models man- 
age to broad ly reproduce th e observed galaxy distribution, 
although the lXu et alJ (l2003f) model predicts the peak to be 
at a higher redshift, of z~ 1. 

The distribution of sources indicate that below a red- 
shift of z~0.5, ~22 per cent of galaxies are early-type (fit 
as ellipticals). Given that our template-fitting classification 
is by nature statistical rather than individual this compares 
favourably with detailed morphological studies of K-selected 
galaxies with deep HUBBLE Advanced Camera for Sur- 
veys (ACS) imaging (Cassata et al. 200^: 20 per cent) and 
the 3.6iim-s elected galaxy population in the ACS field of 
UGC10214 faatziminaoglou et alJl2005at 24 per cent). In 
the IRAC channels, the early-type population extend out to 
a redshift of around unity. Above this redshift the number of 
sources drops away as the sensitivity of the survey to these 
sources is reached. The remainder of the galaxies are spi- 
rals and starbursts. These extend out to redshifts of around 
1.5 before they start to drop below the sensitivities of the 
optical survey. 

The proportion of early-types in the 24u.m population 
drops to ~ 15 per cent of the galaxy population. Tradi- 
tionally, the IR emission from early-types is expected to be 
low since they are old stellar populations with little dust 



1 At points in this paper, rather than present results in all bands, 
we focus only on the 3.6 and 24um results, since they represent 
the two ends of the populations under study. 



or ongoing star formation, leading to some contribution to 
the IRAC bands but negligible contribution to 24u.m (e.g 
IRAS observations of the Shapley-Am es sample found very 
little emission at 60 u.m for ellipticals; Ide Jong et aljll984l : 
Ivan den Bergbll2004l) . Ho wever, IR-lum i nous ellipticals are 
not unknown: for example iLeeuw et alJ (l2004h 's analysis of 
FIR-bright ellipticals sugges ted they cont a ined c old dust in 
their central regions whilst iKrause et "ail (l2003h 's observa- 
tions of one of the brightest YJQ\im sources in the ISOPHOT 
Serendipity Survey showed it to be an elliptical with a pos- 
tulated opaque hidden starburst in its centre. 

The 2A\im population of early-type galaxies in SWIRE 
is therefore an interesting one- it may well turn out that the 
frequency of dusty early-type systems is greater than previ- 
ously thought. As noted in RR05, it is likely that some of the 
elliptical template fits to 24u.m sources will be Arp220-like 
objects which, in the optical, have SEDs similar to elliptical- 
like systems due to the heavy obscuration of young stars. 
Similarly, late-type galaxies with significant dust extinction 
would appear, optically, to be reddened/early-type. A re- 
cent sp ectroscopic study of SWIRE/SPSS galaxies at low 
redshift iDavoodi et alJ 1(2006) has shown a significant popu- 
lation of IR luminous galaxies have red optical colours that 
cannot be explained by optical extinction or contribution 
from AGN. 

The mean redshift of the elliptical, spiral and starburst 
populations in each band can be found in Table Q 

In the right panels of Fig. |S]we plot flux distributions of 
selected IR populations (3.6p.m and 24u:m sources) split into 
redshift ranges. These show a strong evolution from z=0 to 
z~0.8 but this evolution is not seen to continue to higher 
redshifts (as commented by PG05 and lCharv et al.ll2004h . 
These, however, are those Spitzer sources with an optical 
counterpart; the ~30 per cent of sources which are optically 
blank might well be expected to lie at h igher redshift - such 
as in the bimodal redshift distribution of lCharv et al] (I2004T) 
which predicts a second IR bright population centred around 
redshift ~2. 

Type 1 AGN solutions 
Eight per cent of SWIRE sources are identified by photomet- 
ric redshift fitting as optical (Type I) AGN. F05's study of 
faint Chandra X-ray sources in ELAIS Nl found the 24 (im 
population was 10 to 20 per cent AGN (Type 1 and 2) with 
a ratio of between 1:2 and 1:3 between Type 1 and 2. Hence 
~ 5 per cent of the 24[i.m sources were Type 1 AGN. As a 
comparison, of the 18,200 24u.m sources with a good pho- 
tometric redshift solution and 17<r'<23, we find 914 op- 
tical Type 1 AGN fits, or 5 per cent. Although a detailed 
cross-analysis of X-ray-derived AGN fractions to the SED 
template fitting and compact morphology approach of this 
paper is not conducted here, we have made a simple cross- 
match between the photometric redshift catalogue and the 
99 X-ray sources from F05. Of the 50 matches, 16 are classed 
as optical Type 1 AGN from our photometric redshift SED 
fitting whilst F05 class 11 of these 16 as Type 1 AGN, 4 
as Type 2, and 1 as a late- type galaxy. Thus our method 
is successful in identifying optical Type 1 AGN. There will, 
unavoidably, remain some Type 2 AGN 'contamination' in 
our galaxy samples - for example, of the four elliptical-SED 
sources that obtain a match to the F05 Chandra sample, 
the F05 designation was split into early-type (2), galaxy (1) 
and Type 2 AGN (1) classes. More generally, F05 found 
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that a significant fraction (up to 40 per cent) of IR sources 
classified as AGNs in the X-ray do not show any AGN evi- 
dence in the IR/opical SED - consider that all of our 18,200 
24[j.m sources have an IR excess, of which 23 per cent are fit 
in the mid-IR as AGN dust tori (i.e. Type 1 and 2 AGN). 
This implies a Type 1 to 2 ratio of 1:4 or 5 since only 5 
per cent are fit as a Type 1 AGN in the optical. This is 
the canonical fraction one finds from standard X-ray back- 
ground (XRB) models but these tend to be driven by Type 

1 AGN work, which peak at z^2 say, whereas F05 found 
the ratio was reduced due to XRB contribution from Type 

2 sources at moderate to low redshift. Given this we would 
expect the obscured Type 2 to be 'hidden' in the galaxy LF's 
- the reader is cautioned to keep this caveat in mind when 
interpreting the results, particularly at the longer Spitzer 
wavelengths where it is notoriously difficult to disentangle 
AGN/starburst contributions. 

The AGN are seen to be a population that extends 
across the full redshift range, with a slight peak in their 
distribution at a redshift of around two, and a tail to high 
redshift. This slight rise to redshift two, and decline ther e- 
after, is as predicted by the new version of lXu et al.l (|2003)'s 
models when including a new evolutionary model for dusty 
galaxies and the SWIRE limits. This pred i cted AG N distri- 
bution was given in Hatziminaoglo u et al ] J2005bl) . 

Above a redshift of around two the AGN solutions dom- 
inate over the galaxy counts, though this is mainly due to 
SWIRE's ability to detect them to higher redshift. 

3.2 Mid IR SED fitting 

Mid-IR template fitting is carried out in a technique de- 
scribed in RR05. A source that has been fit with an optical 
to near-IR galaxy template, based on its U band to 4.5[J.m 
detections as set out in >I3.1I now has its IR excess calculated 
by subtracting the galaxy model fit prediction from the 4.5 
to 24um data. 

At least two of these bands then need to exhibit an IR 
excess (one of which is required to be 8 or 24[im). This excess 
is t hen characterised by finding the best- fitting out of a cir- 
rus (lEfstathiou fc Rowan- Robinsor]|2003l) . M82 and Arp220 
starbursts jEistathiou^rtal I2000T) or AGN dust torus IR 
template llR^wan^oThrsonT 1^95^ (shown in Fig. More 
information on these templates can be found in RR05. Here, 
it is merely noted that each template is the result of radia- 
tive transfer models and that although for individual galax- 
ies one could obtain better individual fits using radiative 
transfer codes, we here wish to characterise only the broad 
types of IR population. This 'first-order' approach is in fact 
remarkably effective and a blind comparison of predicted 
70\im flux (without making use of the 70[im information) 
derived from the mid-IR SED fitting method adopted here 
to the observed 70fi.m fluxes produces impressive agreement 
for all four mid-IR templates, across almost 2 orders of mag- 
nitud e in 70u.m flux (see fig. 5 of iRowan-Robinson et alJ 
l200fj) . 

Sources are also allowed to be fit by a mixture of 
an M82 starburst and cirrus since it was found that, 
often, both components were required to properly rep- 
resent the IR excess. This mirror s the conclusions of 
Rowan- Robinson fc Cra wfor 1 Jl989ft from fitting mid-IR 
SEDs to IRAS sources, and the findings of RR05. 
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Figure 6. Mid-IR spectral energy distributions: four com- 
ponents - cirrus (split into low-mass [dashed line] and high-mass 
[dotted line] components in the optical region), M82-like star- 
burst [solid line], Arp220-like starburst [dotted line] and AGN 
torus, showing assumed optical/IR ratio at log Lgo=14 (upper 
curve) and 8. Taken from Rowan-Robinson (2001). 



For an IR SED to be fit, then, there needs to be an IR 
excess. If this is not the case then an mid-IR template-based 
K-correction, K/_r, is not calculated for the source and the 
extrapolated mid-IR contribution from the optical/near-IR 
template fit is used instead for in the calculation of z m ax, ir 
(see In the case of a single band IR excess, an M82 SED 
is assumed. Statistics for the IR excess sources are given in 
Tabled 



4 LUMINOSITY FUNCTION PROCEDURE 

Now that the optical/near-IR and mid-IR SEDs of each 
source have been characterised and their photometric red- 
shifts determined, this information can be used to calculate 
LFs, and to include K-corrections derived from the model 
templates. 

Provided a selection function can be defined, a gen- 
eral luminosity function is calculable so long as there is no 
significant population of objects which remain undetected 
at any redshift 2 . Assuming that we have a sufficiently de- 
tailed selection function (see H4.1H and that no such invisible 
population exists - a reasonable assumption given that lo- 
cal galaxies can all be expected to be detected in some of 
the multi-variate waveleng ths under cons ideration - we can 
use the 1/V m ax meth od llSchmidd Il968l and extended by 
lAvni fc Bahcalll Il98fil to determine LFs for each Spitzer 
waveband. 

In this, the LF, $(L), is calculated using the 1/V max 
method where Vmax is the volume corresponding to the 
maximum redshift at which a source could be detected by 
the survey. This maximum redshift may be set by either 
the optical limits (r' < 23), or the IR limits, as determined 
by using K-corrections calculated from the optical/near-IR 

2 There is also a possibility that so few galaxies of a population 
are detected that a selection function cannot be derived for them. 
Some discussion of potential bias is given in i | 1.21 
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Figure 4. X^, ed of photometric redshifts for the different redshift bins considered. Left panel are solutions for galaxies; redshift bins 
of [0 to 0.25] (black solid); [0.25 to 0.5] (green dotted); [0.5 to 1] (cyan dashed); [1 to 1.5] (purple dot-dash); [1.5 to 2] (red triple-dot-dash). 
Right panel arc solutions for AGN; redshift bins are; [0 to 1] (black solid); [1 to 2] (green dotted); [2 to 3] (cyan dashed); [3 to 4] 
(purple dot-dashed). Also shown is the maximum allowed X^ed = ^ cut-off. 




Figure 5. Photometric redshift distributions in selected Spitzer bands. Left panels: N(z) for all sources (black histogram); 
Galaxy template solutions (E, Sab, Sbc, Scd, Sdm and starburst) are overplotted as light blue dot-dashed histogram; ellipticals are 
overplotted as purple dashed histogram; AGN solutions are overplotted as red dash-triple-dot histogram. The mean redshift of each 
population is indicated by a vertical line; light- blue dot-dash line (all galaxies), purple-dashed (ellipticals), red triple-dot-dashed (AGN). 
Right panels: Number counts for sources in several redshift ranges (not corrected for completeness). 
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and mid-IR template fits. We choose to be conservative in 
our choice of IR flux limits (to minimise incompleteness 
and maintain reliability, see i)4.1jl . setting the lower IR flux 
limit below which a source is not considered for a given LF 
calculation to be: 30[jJy (3.6uro); 30^Jy (4.5pm); 50ujy 
(5.8um); 50nJy (8u.m); 500ujy (24(im). 
Thus: 

$(L)dL = J2 77^— 

V max,i 

i 

where "^(L) is the number of objects per Mpc 3 in the (rest- 
frame) luminosity range L to (L+dL) and 



V m 



mm(V r op ticai,i, ViR.i) 



(2) 



So for the case of the rest-frame 3.6u.m LF, Vir,i is ob- 
tained by taking the source in question and redshifting it 
until its observed flux at 3.6(i.m drops below the chosen limit 
of 30/iJy. Vopticai.i corresponds to the maximum redshift out 
to which the same source would be detectable, down to the 
chosen magnitude limit of r'=23 (in order to minimise in- 
completeness and photometric inaccuracy; see Fig.^l. 

The associated rms error, from Poisson statistics, is 
given by 



(3) 



To find Vopticai.i we find the redshift, z m ax,o P tical, at which 
the optical apparent magnitude, m(z), of a source of abso- 
lute magnitude M, when redshifted to a redshift, z, becomes 
fainter then the chosen limiting magnitude of the survey 
(r' < 23), where m(z) is given by; 

m(z) = M + 5log 10 (-^-) + 2.5logi K opt (4) 
Wpc 

Here, Dl is the luminosity distance to redshift z, and K opt 
is the K-correction. K op t also includes the increased effec t 
of IGM absorption at higher redshifts iMadau et aljll993) . 
K opt is defined as: 



A' opt (z) = (1 + 2) 



/ F(\)S{X)d\ 



J F(j±-JS(\)IGM(\,z)d\ 



(5) 



where F(A) is the optical SED template at source, S(A) is 
the combined filter transmission and CCD response curve, 
and IGM(A, z) is the redshift dependent IGM transmission. 

Similarly,for Vma, we find the redshift, z ma x,m at which 
the IR flux in the IR band in question, Sir(z), of a source 
of absolute monochromatic luminosity, Lir, when redshifted 
to an arbitrary redshift, z, becomes fainter than the chosen 
limiting flux of the survey. The K-correction, Kir, is defined 
as in Eqn.|5] with F(A) now the IR SED template at source 
and S(A) is the chosen IR filter transmission. For IR K- 
corrections we do not need to incorporate IGM transmission 
effects for the redshifts considered. 

In order to calculate the LFs correctly, however, we need 
to use a weighting function on the volumes in order to cor- 
rect for the various selection biases that are inherent in the 
catalogue. This is set out in H4.ll 

When calculating the LFs, galaxies with redshifts 
above two are cut since the reliability of their photometric 
redshifts is lower (and the availability of spectroscopic 
redshifts to validate the redshift code is low) and they can 



have a large effect on the bright end of the LF, particularly 
given the low number of sources above this redshift. This 
excludes ~ 0.2% of the sample. Similarly, AGN with 
redshifts above four are cut. Due to photometric variability, 
the reduced accuracy of photometric redshift techniques 
when applied to AGN and the sparser opportunities for 
spectroscopic redshift validation this cut is prudent, exclud- 
ing ~ 5% of the sources fit in the optical by AGN templates. 

We choose to split the LFs into different redshift bins 
in order to probe evolution. For redshift bins the LF cal- 
culation is slightly modified, in that the minimum redshift 
of the redshift bin in question needs to be subtracted from 
the V max ,i in Eqn. [5] The sources are split into galaxies and 
AGN based on their best-fitting optical/near-IR template 
and for these two samples LFs are constructed in the fol- 
lowing redshift bins: 

(i) For galaxies, LFs are constructed in the redshift bins 
[0 to 0.25]; [0.25 to 0.5]; [0.5 to 1]; [1 to 1.5]; [1.5 to 2]. 

(ii) For AGN (smaller number statistics and broader red- 
shift range), LFs are constructed in the redshift bins [0 to 
1]; [1 to 2]; [2 to 3]; [3 to 4]. 



4.1 Incompleteness and selection functions 

The data set used to calculate the luminosity function 
in each IR band is a biased sample. This is due to the 
multivariate flux limits imposed in order for a source to be 
included in the calculations and the requirements that it 
has enough information to be passed to the redshift and IR 
SED fitting codes and obtain a solution (with Xred < 10). 
The general procedure for deriving an overall selection 
function is the same regardless of the IR band in question. 
Here the 24u.m LF will be used as an example: 

The starting point for the 24u.m LF calculation is the 
SWIRE extragalactic catalogue (so stars have been removed 
and 17 <r' < 23) further cut to S(24)>500^Jy (the chosen 
24u.m flux limit from Q. The probabililty that a source of 
a given flux at 24u.m is detected by the survey, ?7detec(£) 
has been determined by Shu pe et al. (2006:in prep.) - for 
IRAC see lSurace et alj |2005). Then, in order to be passed 
to the photometric redshift code it needs to have detections 
in at least four of the bands used to calculate redshift (op- 
tical, 3.6|a.m and 4.5p.m, two of which need to be 3.6\im 
and r'. This is characterised by 7j pz ,i n (5'24). The source then 
needs to obtain a solution with Xred< io> as characterised by 
(Sm). 

Thus the probability that a given 24u.m source has suf- 
ficient information to be passed to the photometric redshift 
code, 7?pz,in(S , 24), is calculated in the following manner: 

We take the input catalogue and construct a source 
count distribution, N a n(S). We then find all those sources 
which are eligible for photometric redshift calculation and 
construct a similar source count distribution, N pZj i n (S). Di- 
viding one by the other gives the selection bias, as a function 
of 24[i.m flux: 



f?pz,in(S) 



(6) 



Wall OS) 

Similarly, we can calculate the selection bias, ?7pz,out(S24), 
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introduced by requiring a source to receive a redshift solu- 
tion with Xred < 10, by constructing the source count dis- 
tribution, N pz ,out(S), of all sources which did obtain a good 
redshift solution, and obtain the resulting selection bias via; 

??pz,out \& ) — 



V pz,out 



(S) 



(7) 



An overall selection bias, 77(5*24), as a function of 24u.m flux 
is then the combination of these functions, namely; 



T](S) = 77dctcc(5')r7pz,in(5')r7p z ,out(S) 



(8) 



The resulting selection functions for each IR band are 
given for IRAC an d MIPS 2A\im in Fig. Following 
ISerieant et alJ i200lt) these selection functions are used to 
weight the contribution of each source to the volume ele- 
ment integral in the LF calculation by modifying the volume 
elements in Eqn.Qto become; 



V' ■ = 

v max. 



f?(z,Si)— dz 



(9) 



It is important to note that this process will not be 
able to fully account for blank-field sources - those sources 
detected by Spitzer but absent from the optical catalogues 
and so having an undetermined redshift distribution. This is 
one reason for the choice of conservative IR flux limits. After 
implementation of the lower flux cut, requiring 17 <r' < 23 
and specifying a photometric redshift solution with Xred < 
10 the number of sources contributing to the LF calculation 
is as follows: 

For 3.6u.m LF calculation - 34,281 sources. 

For 4.5u.m LF calculation - 34,281 sources. 

For 5.8u.m LF calculation - 13,623 sources. 

For 8[im LF calculation - 17,386 sources. 

For 24[im LF calculation - 4,800 sources. 
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4.2 Bias 

LFs are the result of complex selection functions and correc- 
tions for incompleteness required to account for the biases 
and limits introduced by observational selection effects. 

The LFs calculated here are the sum of the LFs of each 
galaxy type. As such, there will be bias if certain galaxy 
types are not visible in a given redshift bin. However, 
even when they are, there will still be a potential bias if 
they are observab l e in different absolu te magnitude ranges 
jLillv et al] Il99& iBrown et all l200lll . There are several 
different LF estimat ors, each with their own advantage s 
and drawbacks (e.g IWillmerl 119971 iTakeuchi etaH l2Q0CfL 
The different bias es each method can introduce to the LF 
were explored by lllbert et al] (120041) who considered the 
dependancy of limiting absolute magnitudes and the intrin- 
sic bias introduced by the differing LFs of differe nt galaxy 
popu lations. The classical 1/V max estimator dSchmidd 
1968), as used in this paper, accounts for the apparent 
magnitude limit of a survey by assuming a homogeneous 
galaxy distribution. It has been s hown that this assump tion 
is not robust to clusters or voids llTakeuchi et al]l2000fl but 
this should be less of an issue for the 6.5 square degree 
field and large redshift range used here. An advantage is 
the simplicity of the method and that it is able to give 
the shape a n d nor malisation of the LF at the same time. 
lllbert et al] i2004h found that the 1/Vmax method acts 
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Figure 7. IRAC and 24|ajn selection functions. The chosen 
flux limit for each band, is shown as a vertical dot-dashed line. 
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Table 1. Statistics from the catalogue. Note these statistics are for the full coverage region, prior to any IR flux cut imposed for 
LF calculations. 





3.6|xm 


4.5|xm 


5. 8um 


Sum 


24um 


Enough bands for z p / lot , N4t, an d 


181,000 


138,000 


40,000 


35,000 


24,000 


After stars removed, N4b an d, e xgal 


140,000 


105,000 


26,000 


25,800 


23,500 


Specify 17<r' <23, N 4bandoxgal-r 


103,000 


81,000 


23,000 


24,000 


20,100 


N4band,cxgal with Z phot solution, N z 


134,000 


99,500 


25,600 


25,500 


22,800 


N 4band,oxgal,r with Zphot solution, N z , r 


102,000 


79,800 


23,200 


24,100 


20,000 


NT*,,, with x d < 10, ^>z,r x 


96,200 


74,700 


20,500 


21,600 


18,200 


^s,TX' galaxies, N aa j and (z) 


88,600; 0.55 


71,000; 0.55 


19,700; 0.50 


20,700; 0.40 


17,300; 0.53 


N z ,rx* ellipticals and (z) 


12,900; 0.38 


12,500; 0.38 


6,390; 0.36 


4,470; 0.30 


1,630; 0.28 


N ZjI . x , [Sab, Sbc, Scd, Sdm] and (z> 


52,600; 0.53 


41,700; 0.54 


11,100; 0.54 


13,400; 0.41 


11,500; 0.52 


N Zj rx* starbursts and (z) 


23,100; 0.67 


16,800; 0.71 


2,260; 0.66 


2,850; 0.51 


4,130; 0.68 


N s ,r x , optical AGN, N agn and (z) 


7,610; 2.78 


3,680; 2.57 


756; 1.99 


880; 2.01 


914; 2.02 


N z rx with IR excess, N excess 


26,500 


26,300 


14,000 


20,600 


18,200 


Neicess fit as cirrus 


4,800 


4,810 


3,440 


4,610 


4,570 


N<, XC ess At as M82 


5,510 


5,500 


3,270 


4,330 


4,230 


N eI „ sa fit as Arp220 


2,670 


2,670 


1,790 


2,590 


2,240 


N eIces5 fit as AGN dust torus 


5,160 


5,150 


3,400 


3,470 


4,290 


N eIces5 single band IR excess 


8,303 


8,190 


2,050 


5,600 


2,840 



to underestimate the LF at the faintest luminosity bins, 
biasing the global LF faint-end slope. This was because 
beyond the absolute magnitude where a certain galaxy 
population is not observable, the method recovers the LF 
slope of the remaining populations. For the 1/V max LF 
calculated here, there is therefore a potential that the 
faint-end of the LFs has been underestimated. In general, 
however, if the selection function as a function of flux and 
redshift is accurate, the 1/Vmax will give an an unbiased 
estimate. 

The second consideration for the LFs calculated here 
is that they are based on photometric redshifts. One ap- 
proach is to assume the photometri c redshift is true for each 
source (as in other studies, such aslOwvn fc Hart,wicklll99fit 
I Sawicki et ail Il997t iMobasher et all I199& iTakeuchi et alJ 
2000). This leads to potential problems with the derived 
LF, which have been minimised to some extent by exclud- 
ing poor xled solutions, only calculating the LFs in the red- 
shift range where the photometric method has been tested 
to a reasonable extent (z of 0-2 for galaxies and z of 0-4 for 
AGN), and by using generous redshift bin widths. Despite 
this, the effect on the LF is a systematic one, tending to in- 
crease the number density at the bright end of the LF since 
the combination of redshift uncertainty and the steep slope 
of the LF at the bright end serves to scatter more sources 
to higher rather than lower luminosity. As a result the LF 
becomes flatter t han the intrinsic LF at the bright end (see 
IChen et al ]|2003l for an in-depth discussion) . In order to ob- 
tain a more quantitative measure of the effect on the LF 
of the photometric redshifts we here carry out Monte Carlo 
simulations where each photometric redshift is assumed to 
be sampled from an underlying Gaussian distribution. This 
will enable the effect of photometric redshift uncertainty to 
be quantified in the LF bins. 

Monte Carlo Analysis 

Each source's photometric redshift is assumed to be sampled 
from an underlying Gaussian distribution with a width, a, 
that is dependent on both the redshift and template type of 
the photometric redshift solution. This dependence is cho- 
sen since the accuracy of the photometric redshift method 
varies with redshift and is poorer for AGN in comparison 



Table 2. Gaussian parameters: The chosen a of the underlying 
Gaussian distribution adopted for the Monte Carlo alteration of 
the photometric redshifts. This a refers to a Gaussian centred on 
the photometric redshift, in log(l+z) space. 



Galaxies 




AGN 


Redshift Interval 


(7 


Redshift Interval a 


0<z p hot < 0.5 


0.05 


0<z p hot < 0.5 0.1 


0.5<z phot < 1.4 


0.1 


0.5<z phot < 2.5 0.2 


Zphot >l-4 


0.15 


z phot > 2.5 0.15 



to galaxies (from the studies in M. II and B04). Consider the 
following: 

(i) Galaxies. At lower redshift, z<0.5, the source will typ- 
ically have good photometry and SED features such as the 
Balmer break are well bracketed by the filters. At higher 
redshift, 0.5<z<1.4, there is still a reasonable bracketing of 
the main SED features by the filters, but the photometry 
is poorer. Hence the accuracy of the photometric redshift 
is reduced, though it has been shown by numerous groups 
that for galaxy samples so far tested the acc uracy of photo- 
metr i c redshift methods is within a z <0.1 dConnollv et alJ 
Il997t iFernandez-Soto et"aH l200ll: RR03; B04; RR05). At 
still higher redshifts, z>1.4, the Balmer break has left the 
z' filter, but the Lyman limit feature will not reach the U 
band until a redshift of 2.5. Coupled with increasingly faint 
sources, the redshift accuracy can be expected to be further 
reduced. 

(ii) AGN. Again, at low redshift, z<0.5, the source will 
be brighter and the photometry well determined (aside from 
photometric variability). At higher redshift, 0.5<z<2.5, the 
photometric accuracy is reduced, mainly due to degeneracies 
between different emission lines. At high redshift, z>2.5, the 
source has become fainter, but the Lyman limit SED feature 
has now entered the filters, improving the redshift identifi- 
cation (i.e. the Lyman dropout technique). 

Based on these considerations, the relevant Gaussian is de- 
fined for each source in Table |2] 

In order to carry out the Monte Carlo analysis each 
source's photometric redshift is replaced by a redshift drawn 
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randomly from the Gaussian distribution centred on the 
original photometric redshift (with the specification that the 
new redshift cannot be less than zero). The usual LF pro- 
cedure is then carried out. This is carried out 100 times, 
for the LF in each Spitzer band, thus producing a series of 
LFs for a given band. The rms spread in these LFs then 
represents the uncertainty introduced by using photometric 
redshifts instead of the true redshifts and this Monte Carlo- 
derived error is found to dominate over the Poisson-derived 
LF error, except in the least populated luminosity bins (the 
highest redshifts and luminosities). The error in the AGN 
LFs will be larger since the photometric redshift technique 
is less accurate for AGN, due to their photometric variability 
over the timescale of the optical survey. 



5 RESULTS 

5.1 Galaxy luminosity functions with redshift 

The galaxy LFs are plotted in panels of Fig. |S| , with errors 
derived from the Monte Carlo treatment. 
The galaxy LF at 3.6u.m and 4.5u.m exhibit the classi- 
cal Schechter functional form JSchechterlll97rj) which re- 
flects the galaxy mass function (e.g. lFont^naet~all l20o3: 
Bund v et al . 2005). That is, the 3.6 and 4.5(i.m LFs are sam- 
pling the stellar mass of the galaxies, since the bulk of emis- 
sion at these wavelengths still arises from the older stellar 
populations. At the longer Spitzer wavelengths the LF is 
then essentially sampling the star formation rate via dust 
emission, and the characteristic break seen at shorter wave- 
lengths is less prominent. The bright upturn (for L>3E10 
i/L^/h^L©) seen in the low redshift, z<0.25, 2A\im galaxy 
LF in Fig. |U can be plausibly attributed to the contribution 
of obscured AGN since, as discussed in >I3.1.3I we expect an 
increasing input from Type 2 AGN to our galaxy LFs at the 
longer Spitzer wavelengths due to the classification criteria 
we use (optical Type 1 AGN or galaxy). 

At the shorter wavelengths these LFs do not demon- 
strate strong evolution where the results are consistent 
with a passive pure luminosity evolution to higher lumi- 
nosity /redshift. The most obvious evolution is seen at the 
longest wavelength (24u-m) where for the first three redshift 
bins this evolution can be seen as a general strong increase 
in luminosity with redshift. The two highest galaxy redshift 
bins (z>l) do not show a strong continuation of this trend. 
Across all five bands we find that the numbers of galaxies 
in the highest redshift bin (1.5<z<2) are not sufficient to 
properly characterise the galaxy LF at this epoch, but the 
results hint at a drop in space-density. 

One might additionally expect that different types of 
galaxy exhibit a different luminosity function due to their 
differing formation and evolution history. We can use our 
broad-brush SED characterisations in order examine the 
galaxy LFs with source type. The photometric redshift 
template-fitting approach allows us to constrain the redshift, 
and to loosely constrain the source type - that is, the char- 
acterisations are valid in a large-scale statistical sense but 
have less meaning for individual sources. The ability of pho- 
tometric redshift codes to reliably recover galaxy type in 
addition to redshift is not well quantified, since the method 
often relies on clear SED features such as the Balmer break 




Figure 8. Galaxy luminosity functions. Top Panel to Bottom 
Panel: LF for galaxies at each wavelength (excluding AGN), split 
into redshift bins of [0 to 0.25] (black); [0.25 to 0.5] (green); [0.5 
to 1] (cyan); [1 to 1.5] (purple); [1.5 to 2] (pink). 
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whic h are found acros s several galaxy types. The study of 
Bolzonclla ct al found, however, that the success rate 

for bursts and ellipticals was high, since they differ markedly 
from more intermediate types. Thus here we choose to split 
the galaxy LFs by optical SED template type, into three 
sub-classes (named so as to avoid morphological inferences); 
Early-type, 'Late-type' (Sab, Sbc, Scd and Sdm), and Star- 
bursts, with this caveat in mind. 

The LFs of these sub-classes are plotted in sub-panels of 
Fig- EI for the 3.6(i.m and 24u.m populations. In each panel, 
the LFs for each sub-class are shown broken down into each 
redshift bin. From this procedure we see that the 3.6[im 
late-type and starburst galaxy populations increase in ac- 
tivity out to z~l. The early- types do not contribute signifi- 
cantly to the evolution, and are consistent with a population 
formed at high redshift and then dimming as they age. The 
evolution in the 24u.m galaxies is also seen to be due to late- 
types and starbursts in the range z~ 0.25 to 1 and at both 
wavelengths the results would support an increase in the fre- 
quency of starbursts (and hence by implication, mergers) in 
the past. 

T 1 e study of lFranceschini et al] ]2006l) combined multi- 
wavelength data in the Chandra Deep Field South, in- 
cluding very deep optical and near-IR photometry, deep 
(S(3.6(j.m) >1 /iJy) IRAC data and complementary opti- 
cal spectroscopy. Importantly, they were able to differenti- 
ate reliably between early- and late-type classes out to high 
redshift and thus confirmed evidence for 'downsizing' occur- 
ring in galaxy formation such that more massive galaxies 
completed their star formation and star accumulation ear- 
lier than less massive ones. Their interpretation was that of 
a morphological shift from the star-forming to the passively 
evolving phase (due to mergers and gas reduction) from high 
to low redshift, something only hinted at by the population 
breakdown shown in Fig. [§] It is of interest to note that 
the 3.6u.m LF for starbursts (bottom-left panel of Fig. |UJ 
shows little evolution at high luminosity and more marked 
evolution at lower luminosity, as would be expected if less 
massive/luminous galaxies remained active to later epochs. 

The depth of the survey does not allow strong conclu- 
sions to be made about the behaviour of early-types above 
z~l. The elliptical template used is of order 10 Gyr in age, 
though stellar synthesis models indicate that this should 
be reasonably similar to galaxie s as young as 4 Gyr af- 
ter an instantaneous bu rst (e.g. iBruzual fc Chariot! Il993t 
iBruzual fc Charlotl l2003h . Galaxies younger than this will 
appear bluer and more like intermediate-type galaxies. Since 
the age of the Universe ranges from ~ 3.1 Gyr to ~ 5.6 Gyr 
in the range 2<z<l we would expect few galaxies to be fit 
well by our elliptical template at these redshifts. 

5.2 Type 1 AGN luminosity functions with 
redshift 

The Type 1 AGN LFs are plotted as a function of redshift 
in each Spitzer band in the Fig. 1101 Since the Type 1 AGN 
are a smaller population and suffer less precise redshifts, 
the resulting LFs are noisier. However, the results are clear 
enough to show positive evolution out to at least the third 
redshift bin (2<z<3), with the strongest evolution in the 
range l<z<2. 

One feature is seen at 24u.m where the LF in the l<z<2 



and 2<z<3 is very similar (as compared to 8u.m where evo- 
lution is seen when moving from l<z<2 to 2<z<3). The 
statistical significance of this would need large-scale spec- 
troscopic analysis for a proper interpretation but is likely to 
be, at least in part, due to the SED shape in these ranges 
- although both 8 and 24u.m are being dominated by hot 
dust emission from an AGN for these redshifts, the 8[im 
emission corresponds to the rising part of the torus just af- 
ter the minimum at Ipum whilst 24u.m traces the peak of the 
torus emission, plus silicate emission. Thus variations in the 
exact SED could lead to differences in the observed LF in 
these bands with redshift, with, for example, the presence 
of the silicate feature in the 24[im band at z~1.3 boosting 
the l<z<2 LF luminosities. 



6 OTHER STUDIES 
6.1 Galaxies 

The IRAC 3.6 and 4.5u-m galaxy luminosity functions 
can be compared to near-IR luminosity functions in the J 
(1.25(i.m) and K 3 (2.2[i.m) bands which, similarly, probe 
the stellar mass in galaxie s . The J and K s luminosity 
functions of IPozzetti et al] |'2003) found no evidence of 
steepening of the faint end of the LFs up to a redshift 
of ~ 1.3 and no evolution of the faint end, whilst finding 
positive luminosity evolution at the bright end of the 
LFs. Their mild evol ution to z^ 1.5 is i n agreement with 
previo us indications <|C owic et al] 119961: ICole et al] liobH : 
ICohen|l2002l: iFeulner et al l 1200311 The more recent 3. 6 urn 
LF of lFranceschini et al] (l2006lL for which half the sources 
had spectroscopic redshifts, extends out to z~1.4 and 
demonstrates a trend for increasing luminosity consistent 
with the dominant stellar populations becoming younger at 
higher redshift. 

Fig. 1111 compares the 3.6[im Franceschini et al. result to 
this work's 3.6[im galaxy LF. The agreement between the 
two works is excellent, and can be considered a validation 
of the photometric redshift and mid-IR fitting approach 
used here. Though not over-plotted for reasons of clarity 
the Pozzetti et al. results, shifted from K s -band to 3.6[im, 
are also in close agreement. 

The 24[im galaxy LF is compare d to the loca l lumi- 
nosity function at 25u.m constructed by|Shurje_j?LalJ dl99S|) 
from IRAS data; to the ISO 15u.m LF o flPozzi et aTT<2004l : 
to the dee p 24u.m LF of PG0 5 and to the spectroscopic 
24u.m LF of lOnvett et~ail fcOOEl) from the SWIRE Lockman 
Validation Field. This comparison is split into three broad 
redshift ranges and LFs have, where necessary, been con- 
verted to 24u.m using an assumed continuum shape of -0.5 
for comparison purposes: 

In Fig. El the 'low redshift', z<0.25, 24^m galaxy LF is 
considered. Shupe at el.'s local 25u.m LF follows the low- 
redshift sample of our 24[im LF well, as does the 15u.m 
LF of Pozzi et al. (median redshift of 0.18) and the low- 
redshift sample from Onyet et al. (z<0.36). The comparison 
to PG05 (z<0.2) demonstrates the complementary nature 
of the datasets used - the deeper PG05 data is able to probe 
the LFs to fainter luminosities whilst the wider SWIRE data 
allows one to go to higher luminosities. In the region of over- 
lap the two results agree well. The bright upturn seen in 
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Figure 9. Selected galaxy luminosity functions - 3.6u.m and 24|J.m. First row; LF for galaxies (excluding AGN), split into redshift bins of [0 to 0.25] 
(black); [0.25 to 0.5] (green); [0.5 to 1] (cyan); [1 to 1.5] (purple); [1.5 to 2] (pink). Second row; LF for ellipticals, split into the same redshift bins. Third 
row; LF for spirals, split into the same redshift bins. Fourth row; LF for starbursts, split into the same redshift bins. 



our z<0.25, 24um galaxy LF for L>3E10 uh v /h~ 2 L can 
be plausibly attributed to the influence of obscured Type 
2 AGN which remain in our 'galaxy' selection and start to 
have an influence at the highest luminosity (particularly at 
these low redshifts where Type 2 AGN are more numerous). 
In Fig. EH the 'intermediate redshift', 0.25<z<l, 24um 
galaxy LFs are compared to the higher-redshift sample of 
Onyett et al. (z>0.36), which exhibits a similar strong evo- 
lution at the bright-end to the intermediate-redshift bins of 
this work. The median redshift of the Onyett et al. sample 
is 0.29, with 38% of the sample in the z>0.36 bin and the 
most distant source at z~2.9. The intermediate redshift LFs 
from PG05 are also plotted. 

PG05's highest redshift results; 0.8<z<1.8, are compared to 
ours; 1.0<z<2.0, in Fig. 1141 Despite the slightly disparate 



redshift bin ranges and luminosities sampled, the two results 
are surprisingly good for these distant populations. 

In summary, then, the SWIRE LFs derived using a large 
photometric dataset have shown good agreement with other 
comparitive works, whilst, crucially, extending the LFs to 
high luminosity at high redshift where the results show con- 
tinued evolution. 

Comparison to GALFORM+GRASIL 

Finally, we compare our 3.6 nm and 24u.m galaxy LFs to 
model predictions from l|Lacev et all 1200.4 Lacey et al. 
2006, in prep.) that combine output fr om a semi-analytic 
galaxy formation model (GALFORM, iBaueh et all 120051 : 
lOole et, all l2000lL with the GRASIL spect rophotometric 
code ijSilva et all Il998t iGranato et all l2000t) . This GAL- 
FORM+GRASIL model firsts calculates the mass assembly 
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Figure 10. Type 1 AGN luminosity functions. Luminosity functions 
for sources fit by IMPZ as Type 1 AGN, split into rcdshift bins of [0 to 1] 
(black); [1 to 2] (green); [2 to 3] (cyan); [3 to 4] (purple). 



and star formation histories of galaxies in the framework of 
the ACDM cosmology, and then computes self-consistent lu- 
minosities for each galaxy, from the far-UV to the sub-mm, 
from which LFs can be readily calculated. 

In Fig. ll5l and Fig. llfcil we compare our 3.6u.m and 24u.m 
galaxy LFs to two models which differ only in the IMF as- 
sumed. In one model, labelled "Kennicutt" in the figures, 
all stars form w ith a solar neighbourhood IMF (we use the 
lKennicuttlll983l parametrization, which is close to Salpeter 
above IMq, but with a flatter slope below this), while in 
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Figure 11. C omparison of 3.6um galaxy LF to similar work. Solid 
circles joined by solid lines is the 3.6iim galaxy LF of this work, split into 
redshift bins of [0 to 0.25] (black); [0.25 to 0.5] (green); [0.5 to 1] (cyan); [1 
to 1.5] (purple); [1.5 to 2] (pink). The results of Franceschini et al. (2006) 
arc shown as stars, for the rcdshift intervals [0.1 to 0.55] (red); [0.55 to 0.9] 
(green); [0.9 to 1.4] (white). 



the other model, labelled "top-heavy" , stars formed in bursts 
triggered by galaxy mergers are assumed to have a top-heavy 
IMF, with a much larger fraction of high-mass stars than in 
the solar neighbourhood. The para meters of the "top-h eavy" 
model are identical to those used in IBaugh et al.l2005l) , who 
found that a top-heavy IMF was essential for their CDM- 
based model to reproduce the observed numbers of faint 
sub-mm galaxies, and have not been adjusted to match the 
data from SPITZER. It can be seen that the IMF models 
diverge at high luminosity and redshift, where the top-heavy 
IMF produces more sources of high luminosity. This is most 
apparent at 24am where more UV photons from the most 
massive stars are re-radiated by dust but even at 3.6um 
there is a clear increase for the top-heavy IMF since more 
high-mass stars evolve though the Asymptotic Giant phase. 
We see that at both 3.6u.m and 24[i.m, our measured LFs are 
consistent with the model LFs based on a top-heavy IMF, 
and lie some distance from the predictions for a standard 
IMF, in particular under-estimating the number of luminous 
galaxies, especially at higher redshift. 

This is a very interesting result. While our measure- 
ments of the LF cannot uniquely determine the form of the 
IMF, it does seem that, when combined with galaxy forma- 
tion models, the observed evolution in the LF at both 3.6u.m 
and 24[i.m favours an IMF in starbursts which is skewed to- 
wards higher masses compared to the solar neighbourhood. 
Such a top-heavy IMF is also suggested by observations 
of chemical abundances in elliptical galaxies and the intra- 
cluster gas in galaxy clusters, which seem to require larger 
yields of heavy elements than are predicted for a stand ard 
IMF iNagashima et all .2005a; Nagashima et al ] l2005bl) . A 
top-heavy IMF would also help significantly in reconcil- 



18 T.S.R Babbedge et al. 



24 micron Goloxies 



10" 



10" 



10" 4 



0<z<0.25 




Pozji et ol.(2004) 
PG05, 0<z<0.2 
Shgpe el ol.(1998) 
Onyell et ol.(2005). 2<0.36 



10" 



10 7 



10" 10 9 



10" 



10" 



y 1 n-3 



10" 



10" 



24 micron Goloxies 



_ Other work: 




PG05, 0.8<J<1.4 
PG05, 1 _•»•«< 1.8 



10 3 



10' 1 10 1 
{» L„/h" 2 L sun ) 



Figure 12. c omparison of 24u.m galaxy LF to similar work: Low 
redshift. The solid coloured lines joining the circles are the LF of this 
work. This is then compared to the 25u.m LF of Shupe et al. (1998); to 
the 15\im LFs of Pozzi et al. (2004) ; to the spectroscopic 24u.m LF of 
Onyett et al. (2005) from the SWIRE Lockman Vaidation Field; and to the 
24u.m LF of PG05. The Onyett et al. spectroscopic LF is not corrected for 
completeness in the faint bins. For the purposes of comparison all points 
have been converted to the same cosmology and shifted to 24u.m (when 
necessary) using an assumed continuum shape of -0.5. The reader is directed 
to the plot for the key. 
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Figure 13. Comparison of 24(j.m galaxy LF to similar work: In- 
termediate redshift. Similar to Fig. now for intermediate redshifts. 



Figure 14. c omparison of 24u.m galaxy LF to similar work: High 
redshi ft. Similar to Fig.tT^l now for high redshifts. 



ing theroretical models of galaxy formation based on CDM 
with results from extragalactic sub-mm surveys. The large 
numbers of distant IR- luminous sources found by SC UBA 
jEales et alJll999l : iBorvs et alJl2003l : IScott et aljl2002D pose 
a significant challenge to CDM-based models, and early ver- 
sions struggled to explain both the observed sub-mm counts, 
and the inferred sta r formation rates in these sources (e.g. 
iGranato et all EoOO') . Later models have therefore invoked a 
variety of solutions to this problem, including a top-heavy 
stellar IMF, as proposed by (iBaueh et alJ2 005'l ; in this way 
the same sub- mm flux can be generated with substantially 
lower star formation rates. Our results lend further obser- 
vational support to this idea which, if true, means that the 
colossal star formation rates (up to several thousand Mq 
per year) currently inferred for distant sub-mm sources may 
need to be revised significantly downwards, pos sibly to levels 
comp arable to those seen in local ULIRGs fe.g. lFarrah et all 



6.2 Quasars 

The evolution of quasars is well established. Most studies 
agree that the space density of quasars above a given lumi- 
nosity increases by a large factor (~ 50) ou t to 2, then 
levels out, and then dec reases to z~5 (e.g. iRichards et all 
l200fJ: ICroom et allErjoi) . The decrease from redshift z~ 2 
to the present day can be interpreted as due to; a reduc- 
tion in the merger rate; a decrease in the cold gas available 
for fueling; and an increase in the timescale of gas accretion 
jKauffmann fc Haehnelt]|2000l) . 

In this work we see a strong increase in the space den- 
sity at a given luminosity from the 0<z<l to the l<z<2 
bin for all bands and further increase in the 2<z<3 bin. At 
higher redshift the number statistics become poorer but are 
consistent with a constant /rising space density. The recent 
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Figure 15. Comparison of 3.6|im galaxy LF to model. The shaded 
regions arc the observed galaxy LF at 0<z<0.25 (orange) and 0.5<z<l 
(lilac), with the dotted line marking the LF and the region representing 
the ltr errors. Overplottcd arc the models of Laccy ct al. for the rcdshifts 
(orange) and 1 (blue), split into Kennicutt (1983) and Top-Heavy IMF 
versions as striped and solid lines respectively. 
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Figure 16. c omparison of 24um galaxy LF to model. The shaded 
regions arc the observed galaxy LF at 0<z<0.25 (orange) and 0.5<z<l 
(lilac), with the dotted line marking the LF and the region representing 
the la errors. Overplottcd arc the models of Lacey et al. for the redshifts 
(orange) and 1 (blue), split into Kennicutt (1983) and Top-Heavy IMF 
versions as striped and solid lines respectively. 
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Figure 17. Comparison of 24|xm Type 1 AGN LF to similar works. 

Comparison to the 15nm ELAIS AGN LF of Matuto ct al. (2006). For 
the purposes of comparison all points have been converted to the same 
cosmology and shifted to 24fxm using an assumed continuum shape of -0.5. 



8um type 1 quasar evolution of lBrown et all ll200d) . derived 
from S(24um)>lmJy sources with l<z<5, found the peak 
quasar space density to be at z=2.6±0.3 whilst here the 
8um luminosity evolution shown in Fig. I19l rises strongly to 
2<z<3 and then flattens off. 

Hence our broad-brush photometric redshift approach 
finds a similar result to other quasar studies whilst extend- 
ing the Type 1 AGN LFs to the Spitzer wavebands. How- 
ever, due to the depth of the SWIRE survey we are unable 
to fully observe the onset of strong decline at the highest 
redshifts. From optical studies this strong decline in Type 1 
AGN activity occurs above redshifts of ~3-4, which is our 
highest redshift bin, and has be en seen in numerous q uasar 
studies since the early '80 ' s ,e.glOsmedll982UWarren et al 



11994 ISchmidt et aiJll995t iKennefick et alJll995t iFan et aT 



2001). This dramatic drop can be interpreted as marking the 
epoch of peak formation of quasars, at around ten per cent 
of the current age of the Universe. As such, it has strong 
implications fo r models of black hole formation, fuelling and 
evolution (e.g. [Haehnefifc RgesJ^93j: [Granat o et al.ll200ll : 
lHatziminaoglou et alJl2003t iBorvs et alJl200aL 

In Fig.EUthe 24 um T ype 1 AGN LF i s com pared to the 
15um ELAIS AGN LF of iMatute et all J2006h (consistent 
with the earlier LF of lMatute et all2002l, . The Matute et al. 
LF is split into a z=0.1 population and a z=1.2 population, 
which trace a similar behaviour to the evolution with red- 
shift found in this work. The broad redshift bins of this work 
make a direct comparison problematic since wider redshift 
ranges have the effect of flattening the resulting LF, how- 
ever the results are consistent. For closer comparison, the LF 
from this work has been re-binned into narrower redshift in- 
tervals, though the resulting low numbers of 24um Type 1 
AGN in the bins does result in upper limits, particularly at 
higher luminosity and for the lowest redshift bin (z<0.4). 
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7 LUMINOSITY EVOLUTION, ENERGY AND 
SFR DENSITY 

7.1 Schecter and double power-law fits 

The LFs calculated here demonstrate the evolution of both 
galaxies and Type 1 AGN with redshift. We now model these 
LFs using Schecter Functions dSchechte and double 

power-law fits. The Schecter parameterisation is often ap- 
plied to UV/optical galaxy LFs, whilst double power-law 
fits have been f ound to be a good fit to the local infrared 
population (e.g. lLawrence et al ] ll986tlSaunders et aljll990t 
iTakeuchi et al]l2003h . This can be broadly explained since 
star formation acts as a power-law in the luminosity func- 
tion and stellar mass acts as an exponential drop-off - a 
luminosity function can then be seen as a convolution of 
the two regimes, weighted by the relative dominance of each 
(determined by the rest-frame wavelength). 

The T ype 1 AGN LF i s often modelled b y a d ouble 
power- law feovle et alJl988h . For example. iPeil il995l) used 
a double power-law model where the bright and faint ends 
were independent of redshift, with the break occurring at a 
luminosity that was dependent on redshift. 

The Schecter parameterisation is as follows: 



®(L)dL/L* =$*( — ) e 



L 



-L/L* 



dL/L", 



whilst the double power-law is described by: 

l-Q 



$(L)dL/L* = $* 



and 



${L)dL/L* = $* 



L 



dL/L*,L< L* 



1-0 



dL/L* ,L > V 



(10) 



(11) 



(12) 



Our approach to fitting the LFs is determined by the 
wavelength under consideration and the extent to which the 
data can break parameter degeneracies. Thus, we choose to 
split the analyses in the following manner: 



7.2 5.8 and 8um galaxies, and Type 1 AGN LFs - 
luminosity evolution 

Double power-law functions are fit to the 5.8 and 8\im 
galaxy LFs and to the Type 1 AGN LFs (note that Schecter 
fits were attempted for the 5.8\im galaxy LF but these did 
not give good results, most obviously at the bright end). 
However, the luminosity ranges sampled by these SWIRE 
LFs do not allow us to constrain the faint-end well (and 
hence the faint-end slope and the break where the bright-end 
slope takes over). Thus, noting that the bright-end slopes are 
well described by the same slope across the redshift bins, we 
choose to explore the evolution of L* only, by fixing a, j3, and 
$*. We set the faint-end slope to a = 1.3 (cf. PG 05's slope 
of 1—1.3 for their 2A\im sample: Izheng et alJl2006ts value of 
1.2L0.3 from analysis of the 24um LFs of iLe Floc'h et alJ 
2005 and the local val ue for &Q\im IRAS gala xies of 1.23 from 
Takeuchi et"aH 120031 and 1.09 from ISaunders et aljll99cl) . 
The bright-end slope for the 5.8 and 8\im galaxy LFs are 
set by the slope in the lowest redshift bin (0<z<0.25) where 
it is best-constrained, giving /3=2.5 and 2.6 respectively. For 
the Type 1 AGN LFs a shallower slope of (3=2.2 is able to 
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Figure 18. L* (L Q h~ 2 ) for 5.8 and 8 u.m galaxies. 5.8 and 8|jm pa- 
rameters result from PLE power-law fits. The overlaid dotted lines, renor- 
malised to the 5.8iim value in the first redshift bin, show the trend if 
L*<1>* fx (1 -f- z)^, for various 7. The best-fit 7 arc 7=1.6_q'q at 5.8\im 
and 7=1.2_„'g at 8|xm (la errors in the parameter space). 



give a good fit across all the bands. The density param- 
eter $* is set to 5E-2 h 3 dex _1 Mpc" 3 for the 5.8 and Sum 
galaxy LFs and to 5E-5 for the Type 1 AGN. It is important 
to note that although we have chosen reasonable values, for 
the purposes of exploring pure luminosity evolution (PLE) 
the exact choice of the fixed a and <£»* values only effects the 
overall scaling of L*, not its evolution as seen by SWIRE. 
The resulting PLE trends are then characterised by 



(i+zT 



(13) 



Fig. ll8l shows the evolution of L* for the 5.8 and 8\im galax- 
ies, whilst Fig. H9l shows the evolution for the Type 1 AGN. 
For the two galaxy LFs L* is seen to increase out to the 
0.5<z<l resdhift bin, and then flattens off at higher red- 
shift. The best-fit 7 is found to be 7=1. 6lg 6 at 5.8um and 



T.2^0 5 at 8\xrn, with the majority of the evolution occur- 



ring since z~l. In contrast the Type 1 AGN show positive 
evolution in all bands, out to the highest redshifts (3<z<4). 
There is a noticeably higher relative luminosity across red- 
shift for A >5vim in comparison to A <5\im. The best-fit 
7's are 7=1.31^ at 3.6um, -y=1.0±^;^ at 4.5nm, 7=3. 8± ;| 
at 5.8nm, 7=3.4±° j at 8um and 7=3.01^ at 24^im. The 
strong evolution a t longer wave l ength s is consistent with 
other studies - e.g. iMatute et alJ <|2006T) found a pure lumi- 
nosity evolution model with Loc(l+z)~ 2,9 was a good fit to 
their 15^m ELAIS AGN LF. 



7.3 3.6 and 4.5[im galaxies - energy density 

We choose to fit Schecter functions to the IRAC 3.6 and 
4.5u.m galaxy LFs. As the data are sufficient to constrain the 
parameter space we allow a, $*, and L* to vary, except for 
the l<z<1.5 bin where we use the a value fit to the 0.5<z<l 
slope due to the reduced data at this high redshift. We are 
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Figure 19. L* (L Q h -2 ) for Type 1 AGN. L* values result from PLE 
power-law fits. The overlaid dotted lines, rcnormaliscd to the 3.6\im value 
in the first rcdshift bin, show the trend if L* oc (1+z)^, for various 7. The 



best-fit 7 arc, respectively, 7 — 1.3^q'^, 
(1(7 errors in the parameter space). 
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then able to integrate the resulting Schecter functions over 
luminosity in order to determine the total energy density, Q, 
provided by 3.6 and 4. 5 urn galaxies at each redshift. This 
integral is plotted in Fig. I2UI showing the contribution from 
the observed data, and from the Schecter-fits extrapolated 
to luminosities not sampled by the observed LFs. Consider 
that an old passively evolving population is characterised by 
PLE, dimming more quickly than a no-evolution model in 
the K band, say (e.g. lEllis fc Jonesll2004t Ide Propris et alJ 
1999). This scenario is essentially what is found for IRAC 
3.6 and 4.5um, which probe the older stellar populations (in- 
deed a PLE Schecter fit describes the change with redshift 
almost as well as one in which both L* and <£>* are allowed 
to vary). Fig 1201 shows that for both bands, the energy den- 
sity gradually increases out to z~0.5-l but flattens, or even 
declines, at higher redshift. 



7.4 24um galaxies - energy and SFR density 

For the 2A\im galaxies we choose to fit a double power 
law. Again, the data is not able to constrain the faint-end 
slope and the turnover and we fix the faint-end slope to 
q=1.3 and $*=5E-2 hMex^Mpc -3 (as before in 
However, as shown in H6.ll the agreement with the PG05 
LF (which is at fainter luminosities) is good and this choice 
of a and $* gives double power-law fits to our data that 
also match well with the PG05 data at the faint-end and 
turnover location. This gives us confidence in extending 
the analysis to calculate the 2A\im energy density, Q.24, by 
summing the contribution of our data and the integral of 
the double power-law outside of the luminosity range in 
which we have data (as in >I7.31 . For the 24um analysis we 
now go further and use £224 to obtain an estimate of the 
total (8-1000u.m) IR energy density, £Iir, and the resulting 
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Figure 20. Integrated energy density f2 (LgMpc 3 ) for 3.6 and 
4.5|4m galaxies. S^3.6 and SI4 5 (black circles) result from the combination 
of the observed data (red triangles) and the integrated Schcctcr-fits outside 
of the luminosity range where there was data (green diamonds). 



implied cosmic-SFR density: 

Total IR luminosity 
The correlation between rest-frame 12-15um and total 
IR luminosity for local IR-lumin ous galaxies has been 
demonstrated numerous times fe.g. IPapovich fc BelJ [20021 : 
ICharv fc Elbazll200l|) to be good, with a scatter of around 
0.15 dex. iBell et alJ <2005ft calculated a conversion from 
observed-frame 24am to total IR l uminosity using t he full 
range of model template spectra of iDale et al.l l|2001), with 
a mean correction factor of, to within 20 per cent, 



L IR ~ WuL v (2i\im). 



(14) 



The mean template spread around this mean was larger, 
0.3 dex. Similar res ults have been found by the study of 
200 1|) for ISO and IRAS galaxies, and by 
2005) (see their figure 8). Thus, we convert 
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our Q24 to SIir using Ecm. 1141 which should be accurate to 
a factor 2 or so jBell et al.H2005l) . 

Star formation rate density 
We use Qir to o btain the SFR density using the relation of 
lKennicutdJl998l) : 



SFR(M yr" 



1.71xlO" 10 L/fl(L e 



(15) 
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Figure 21. Integrated IR energy density fl IR and SFR den- 
sity from 24(im galaxies. (black circles) results from the combi- 
nation of the observed data (red triangles) and the integrated two-power 
fits outside of the luminosity range where there was data (green dia- 
monds) to obtain fl24- This is then converted to O/jR (LqMpc - ^) and 
SFR (M0yr -1 Mpc~ 3 ) as detailed in the text. The overlaid black-dotted 
lines, rcnormaliscd to the value in the first redshift bin, show the best-fit 
trend and la spread if Qir oc (1 + z)^, with best-fit r y— For com- 
parison the Qjr results from PG05 and Lc Floc'h et al. (2005) arc also 
plotted. 



8 DISCUSSION AND CONCLUSIONS 

This paper has extended the photometric redshift method 
to incorporate near-IR data (sub 5[im) in the solution, im- 
proving both the accuracy and reliability of the solutions. 
Application of the extended photometric redshift technique 
to a 6.5 square degree area in the ELAIS Nl field of the 
SWIRE survey has shown that the majority of the opti- 
cal SWIRE population can be characterised as star-forming 
galaxies with a mean redshift of (z) ~ 0.6, with a tail ex- 
tending out to redshifts of z~ 1.5. Early-type galaxies (fit as 
ellipticals) make up around a fifth of the population, with a 
mean redshift of (z) ~ 0.4 and a tail to a redshift of around 
unity. The proportion of 24u.m sources that are early- type is 
seen to be lower than that of the shorter wavelength Spitzer 
sources, however the results suggests that the frequency of 
dusty early-types is higher than pre-Spitzer studies have 
found. Some of these sources could in fact be heavily ob- 
scured Arp220-like objects. In addition to the galaxies, there 
is a population of optical Type 1 AGN (~ 8% of the total) 
that extends across the full redshift range, with indication 
of a slight rise in their distribution to redshift two to three, 
and a decline thereafter. 

Using the photometric redshift solutions, along with 
selection functions designed to account for the bias in the 
multi- wavelength catalogue, allowed luminosity functions to 
be calculated in each Spitzer band, for both galaxies and 
quasars. The relevant K-corrections in the optical and IR 
were derived using the optical/near-IR template solutions 
of the redshift fitting procedure, and separate mid-IR 
template fitting solutions applied to the longer Spitzer 
wavebands. Importantly, this SWIRE sample has enabled 
a joint analysis of galaxy and Type 1 AGN LFs using the 
same sample, selection and method for both populations. 



As discussed in PG05, this conversion results in further un- 
certainty - the SFRs are subject to photometric redshift un- 
certainties, differences between different model templates, 
and the validity of extrapolating low-redshift relationships 
such as that of Eqn. 114) to higher redshift. Thus the result- 
ing SFRs can be expected to be subject to an additional 
systematic error of up to a factor 2-3. 

The cosmic IR energy density, £Iir (left axis), and SFR 
density (right axis) are plotted in Fig. 1211 showing the con- 
tribution from the observed data, and from the two-power 
fits extrapolated to luminosities not sampled by the observed 
24u.m LF. For clarity the additional systematic uncertain- 
ties resulting from using Eqn. I14| and 1151 for conversion are 
not included in the plot. In contrast to IRAC, strong positive 
evolution, characterised by Qjh oc (1+z) 7 with 7=4.5^06; is 
seen for the 2A\im galaxies, with the majority of this evolu- 
tion occurring since z~l. For comparison, Qir results from 
PG05 and Le Floc'h et al. (2005) are also plotted, show- 
ing the same strong evolution. The evolution in star forma- 
tion rate density PG05 measured was oc 

(1+z) 4.0±0.2 Qut tQ 

a redshift of 0.8, consistent with the evolution found in this 
work, and with t he ev olution of 7=3.91^* to z~l, found by 
iLe Floc'h et alJ i2005tl for their 24|xm-selected sample. 

Note that parameter fit values with redshift resulting 
from *J7]are given in Appendix lAl 



The 3.6 and 4.5u.m galaxy LFs, where the LF is probing 
the stellar mass, showed overall evolution consistent with 
the passive ageing of the dominant stellar population up to 
z~1.5. The comoving luminosity density was found to evolve 
passively, gradually increasing out to z~0.5— 1 but flattening, 
or even declining, at higher redshift. Analysis with templatel 
type does, however, offer tantilising support for 'downsizing' 
with starbursts showing little evolution at high luminosity 
and more marked evolution at lower luminosity, as would be 
expected if less massive/luminous galaxies remained active 
to later epochs. 

More pronounced evolution was seen at the longest 
wavelength; the 24u.m galaxy LF, more indicative of star 
formation rate and/or AGN activity, undergoes strong pos- 
itive evolution, with the derived IR energy density and SFR 
density oc (1+z) 7 with 7=4.5lg'g an d the majority of this 
evolution occurring since z~l. This strong evolution was in 
agreement with other equivalent work, however the SWIRE 
LFs were uniquely able to show that this evolution contin- 
ues out to high luminosities at high redshift (where many 
models start to differ). 

In addition to other observational work, comparisons 
to the predictions of a combined semi-analytic and spec- 
trophotometric code were made. While we cannot provide 
meaningful constraints on the exact form of the IMF 
preferred by our LFs, an IMF skewed towards higher mass 
star formation in bursts compared to locally was preferred. 
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As a result the currently inferred massive star formation 
rates in distant sub-mm sources may require substantial 
downwards revision. 

Although the results for the optical Type 1 AGN LFs 
are subject to greater uncertainty due to reduced photo- 
metric redshift accuracy and smaller population size, the 
analysis was able to demonstrate clear positive luminosity 
evolution in all bands, out to the highest redshifts (3<z<4). 
Modelling as L* cx(l+z) 7 gave 7 ~ 1 — 1.3 for 3.6 and 4.5|im 
and stronger evolution at the longer wavelengths (5.8, 8 and 
2A\im), of 7=3.8tS;3, 3.4±°'*and 3.0±g;| respectively. 
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APPENDIX A: LF FUNCTIONAL FIT 
PARAMETERS 



Table Al. Schecter fits, for 3.6 and 4.5pm galaxy LFs. 



Sample 


a 


L*(L©n- 2 ) 


**(h 3 Mpc- 3 dex- 1 ) Redshift 


3.6 pm 


-0.9 


3.6xl0~ 3 


5.6|Jxl0 9 0-0.25 




-1.0 


4.5xl0~ 3 


7.9* : fxl0 9 0.25-0.5 




-0.9 


7.0xl0~ 3 


6.8f : l,xl0 9 0.5-1.0 




-0.9 


4.0xl0~ 3 


7.l£; 9 xl0 9 1.0-1.5 


4.5 pm 


-0.9 


4.5xl0~ 3 


2.8 3 ,- 2 xl0 9 0-0.25 




-1.0 


6.7xl0- 3 


3.2 3: |xl0 9 0.25-0.5 




-1.1 


7.7xl0~ 3 


3.2 3: |xl0 9 0.5-1.0 




-1.1 


6.9xl0~ 3 


3.2 3 ;|xl0 9 1.0-1.5 


Table A2. 


L* 


fit parameters 


, for 5.8 and 8pm galaxy LFs. 



Double power-law, fixed parameters are: a = 1.3, /3=2.5 (5.8pm) 
and 2.6 (8pm), $*=5E-2 f^dcx^Mpc- 3 . Best-fit 7 is given for 
L* oc (l+z)T. 



Sample L*(L h~ 2 ) Redshift 



5.8pm 



2.5*;°xl0 8 
4.5^'fxlO 8 
7.1 4 1 6 2 xl0 8 
5.6|'|xl0 8 



0-0.25 
0.25-0.5 
0.5-1.0 
1.0-1.5 



1.6 



+0.7 
0.6 



8.0pm 5.0§'^xl0 8 

8.9i 4 3 lxl ° 8 
8.9i^i x l0 8 

1.0j;fxl0 9 



0-0.25 
0.25-0.5 
0.5-1.0 
1.0-1.5 



1.2 



+0.4 

0.5 



Table A3. L* fit parameters, for Type 1 AGN LFs. Dou- 
ble power-law, fixed parameters are: a = 1.3, (3=2.2, <3?*=5E-5 
h 3 dex _1 Mpc _3 . Best-fit 7 is given for L* oc (1+z) 7 . 



Sample L*(L0h 



Redshift 



3.6 pm 


2.8 3 .- 2 xl0 9 
5.6j?' 3 ,xl0 9 

l.lJjxlO 10 


0- 1 

1- 2 

2- 3 

3- 4 


1 o+O.l 
^-O.l 


4.5 pm 


2.5 3 .; 2 ,xl0 9 
4.5|'° 1 xl0 9 
6.3^xl0 9 
7.1 8 ; 9 xl0 9 


0- 1 

1- 2 

2- 3 

3- 4 




5.8 pm 


2.24;0 X 10 1() 


0-1 


0+0.4 
^•"-o.s 




1.6 2 |xlO n 
4.07 :9 xl0 n 
2.0 3 ; 2 xl0 12 


1-2 






2- 3 

3- 4 




8.0pm 


2.5f;|xl0 n 
2.0 3 ' 2 xlO n 
5.6i°2 xl0 n 
7.9| : 6xl0 n 


0- 1 

1- 2 

2- 3 

3- 4 


4 +o.4 


24 pm 


1.4 2 - 2 xl0 1() 
LlJ^xlO 11 
1.4jfxlO n 
3.6|;gxlO n 


0- 1 

1- 2 

2- 3 

3- 4 


n+0.5 
d.U_ Q 4 
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Table A4. Double power-law functional fit parameters, for 

the 24um galaxy LF with a = 1.3 and <I>*=5E-2 Ir^dcx^Mpc- 3 . 
Best-fit 7 is given for L* <x (1+z) 7 . 



Sample /3 L*(L h~ 2 ) Redshift 7 



24 um 


-2.5 


3.6f|xl0 8 


0-0.25 


1 6 +0 ' 7 




-3.0 


1.6i|xl0 9 


0.25-0.5 






-3.0 


3.2f : fxl0 9 


0.5-1.0 






-3.0 


4.5|;gxl0 9 


1.0-1.5 





